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ABSTRACT 



This dissertation presents the first theoretical investigation of the Lamb shift in a 
light-front hamiltonian approach: the dominant part of the sphtting between the 2Si 
and 2Pi energy levels in hydrogen is calculated. Also presented for the first time is 
an analytic calculation in a light-front hamiltonian approach of the singlet-triplet spin 
splitting in the ground state of positronium through order a'^. 

We study the QED bound-state problem in a light-front hamiltonian approach. 
We start from a canonical QED Hamiltonian, and set up a general formalism for 
deriving the effective Hamiltonian Hx to some prescribed order in a (with a <^ 1). 
Hx is renormalized by requiring it to satisfy coupling coherence. Then we use bound- 
state perturbation theory (BSPT) to compute the low-lying spectrum of interest in a 
consistent set of approximations to some prescribed order in a and a ln(l/Q;). The 
general formulas are applied explicitly to the positronium and hydrogen systems. 
Renormalization is carried out through order e^, and a nonrelativistic limit of the 
theory is taken: |p| <^ m , where |p| is a typical electron momentum and m is the 
electron mass. Also, in order to derive the results in the few-body sector of interest — 
|ee) for positronium and \ep) for hydrogen — we require our final cutoff to satisfy 
ma^ <^ A <^ ma. This upper bound is the dominant energy of emitted photons, and 
this lower bound ensures that we do not remove the non-perturbative energy scale of 
interest with our renormalization group transformation that we run perturbatively. 

ii 



To my family 



iii 



ACKNOWLEDGMENTS 



A special thanks to my adviser Robert Perry for his guidance, wisdom (especially 
regarding renormalization) and friendship. 

I would like to thank Brent H. Allen for helping me to learn what quantum field 
theory is all about. 

I would like to thank Stan Glazek for useful discussions and for collaboration on 
our positronium paper. 

A special thanks to Koji Harada for introducing the idea behind the Lamb shift 
calculation, and also for stimulating conversations over the past few years. 

Thanks to my fellow nucs Martina Brisudova, Werner Koepf, Lisa Kurth, David 
Makovoz, John Rusnak, Sergio Szpigel and Trey White for useful discussions, and 
Bunny Clark and Dick Purnstahl for their guidance. 

There are many coneheads with whom I have had profitable discussions. I would 
like to mention Edsel Ammons, Stan Brodsky, Matthias Burkardt, John Killer, Gerry 
Miller, Steve Pinsky, Wayne Polyzou, Dave Robertson (who really helped clarify 
renormalization to me over and over and ...) and Ken Wilson. 

I would like to thank Tom Humanic, Furrukh Khan, Steve Pinsky and Junko 
Shigemitsu for serving on my committee. 

And last but certainly not least, thanks to my wife Saeko for all her patience and 
love, and my parents and siblings for their continued love and support over the years. 

iv 



VITA 



June 28, 1968 Born - Powell, Wyoming 

1987-1989 A.S., Western Nebraska Community 

College, Scottsbluff, Nebraska 

1989-1992 B.S., Department of Physics, Colorado 

School of Mines, Golden, Colorado 

1992- 1993 Graduate Fellow, The Ohio State Uni- 

versity, Columbus, Ohio 

1993- present Research Assistant, The Ohio State 

University, Columbus, Ohio 



PUBLICATIONS 

Research Publications 

Billy D. Jones, "Singlet-triplet splitting of positronium in hght-front QED," in 
25th Coral Gables Conference on High Energy Physics and Cosmology, Proceedings, 
Miami Beach, Florida, edited by Behram N. Kursunoglu, Stephen Mintz and Arnold 
Perlmutter, (Plenum Press, New York, 1997). 

Billy D. Jones, Robert J. Perry and Stanislaw D. Glazek, "Analytic treatment of 
positronium spin splittings in light-front QED," Physical Review D 55, 6561 (1997). 

Billy D. Jones and Robert J. Perry, "Lamb shift in a light-front hamiltonian ap- 
proach," Physical Review D 55, 7715 (1997). 



V 



FIELDS OF STUDY 

Major Field: Physics 

Studies in light-front field theory: Professor Robert J. Perry 



vi 



TABLE OF CONTENTS 



Abstract 

Dedication .... 
Acknowledgments 

Vita 

List of Figures . 
Chapters: 



1. Introduction 



2. Light-front field theory and canonical QED Hamiltonian 



2.1 Light-front held theory . . . 

2.2 Canonical QED Hamiltonian 



T. Light-front hamiltonian bound-state problem 





3.1 


Ceneral computational strategy 




3.2 


Step one: Derivation of Hx with G-W transformation 




3.3 


Step two: Diagonalization of Hx 




3.4 


Wegner transformation 




3.5 


Bloch transformation 




3.6 


Renormalization issues in light-front field theory 






3.6.1 Relevant, marginal and irrelevant terminology (on the light- 


front) illustrated through simple examples 




3.6.2 Simple example to elucidate coupling coherence 



vii 



4. Positronium's ground state spin splitting 



4.1 Derivation of Hx tiirougii second order] 

4.1.1 Kenormalization issues in positronium 

4.1.2 H\ tlirougli order e"': exchange and anniJiilation channels| . . 
O Diagonalization ot Hx 

4.2.1 ?-^o, a coordinate change and its exact spectrum 

4.2.2 BSPT, renormaUzation and a Umit 

4.2.3 Calculation of 5AT^ 

4.3 Singlet-triplet splitting: A simpler method 



5. Lamb shift of hydrogen 



5.1 Theoretical framework 

5.2 Lowest order Schrodinger equation . . , 

5.3 Lamb shift calculatimi 

5.3.1 Low photon-energy contribution 

573.2 High photon-energy contribution 

5.3.3 Total contribution 

5.4 Lamb shift summary and discussion . , 



5. Summary and Discussion 



Appendices: 



A. Summary of light-front conventions 

W. Matrix elements of the canonical Hamiltonian used in the dissertation . . 



U. 4th order Bloch and 3rd order G-W 



D. Hyperspherical harmonics and I'bck coordinate change 



D. -vs- Bj^ in positronium as an illustration 



viii 



h\ Averaging over directions! 141 



Bibliographyl 145 



ix 



LIST OF FIGURES 



Figure Page 

2.1 Feynman diagram for the tree-level annihilation amplitude in 0'^ theory] 8 



O (a) 'i'ime-ordered diagram for the tree-level annihilation amplitude in 



theory, (b) 'I'ime-ordered diagram containing vacuum mixings, ft 



vanishes in light-front held theory. 



4.1 Momenta and spin label conventions for the second-order effective elec- 



tron self-energy interaction.] 55 





4.2 This illustrates the spin and momenta label conventions used for this 


Dositronium calculation 












4.3 dM^ is the part of the ground state spin splittings from first-order 



structure constant. The state labels 1, 2, 3 and 4 are explained in 
Eq. ( [4.98|) . The two upper most levels should coincide in a rotationally 
invariant theory; and after including second-order BSPT, they do. . . 77 





4.4 I'he combined ground state spin splitting from hrst- and second-order 


BSP'i' in positronium through order a'^ is illustrated using the same no- 


tation as in i'lgure 4.^^. 5M'^ is given by Eq. (4.99 


) and is calculated in 



Subsection [4.2.3| . The final combined result (on the right) corresponds 

to a rotationally invariant theory 79 



5.1 I'he ehective interactions that add to give the (Joulomb potential. ''H'^ 



implies that the photon energy is greater than A. "L" implies that the 



electron kinetic energy is less than A. We choose ma^ ^ A ^ ma. 



these "H" and "L" constraints can thus be removed to leading order. 



and we are left with — 4e^/q^, the Coulomb potential.] 100 



X 



5.2 Low photon-energy contribution to tlie Lamb sliift. Diagram LI repre- 


sents the shift arising from treating photon emission below the cutofi 


A in second-order BSP'L, where the intermediate electron-proton are 


bound by the Coulomb potential 


(scattering states must be included 


too of course). Diagram L2 is an 


effective self-energy interaction plus 


counterterm (shown in i'lgure 5.3 


) treated in first-order BSPT. P{2,1) 



is the average excitation energy of the n = 2 levels; see Eqs. ( |5.113| ) 

and ( p.ll4[ ) and the discussion above them for details 107 



5.3 


The sum of an effective self-energy interaction (arising from the re- 


moval of photon emission above the cutoff A) and a counterterm. The 


counterterm is fixed by coupling coherence as in Subsection 4.1.1 


. The 




result is the interaction in diagram L2 of Figure |5.2| 




5.4 


High photon-energy contribution to the Lamb shift. These are 


third 


and fourth order effective interactions treated in first-order BSPT. 




These effective interactions arise from the removal of photon 


emis- 



sion above the cutoff A. '6' is an arbitrary scale, required to satisfy 



ma <^ h <^ that was introduced to simplify the calculation. Note 



the 6-independence of the result. The total contribution is a sum ol 



Figure 5.2 and Figure 5.4. Note the A-independence of the combined 



result 115 



5.5 The experimental n = 2 hydrogen spectrum: fine structure, Lamb shift 


and hyperhne structure. F = L + Se + Sp. 





123 



xi 



CHAPTER 1 



INTRODUCTION 



There is much effort being put into solving for the hadronic spectrum from first 
principles of Quantum Chromodynamics (QCD) in 3+1-dimensions using a light-front 
hamiltonian approach IQ. However, low-energy QCD is challenging, and a realistic 
analytic calculation may be impossible. There is a need for exact analytic calculations 
that test and illustrate the approach. The core of this dissertation provides two 
examples of this in Quantum Electrodynamics (QED): analytic calculations of (i) 
positronium's ground state spin splitting through order a"^ and (ii) the dominant 
part of the Lamb shift in hydrogen through order ln(l/a). These calculations are 
based on previously pubhshed work, and respectively. The specific framework 
of calculation, a hamiltonian light-front renormalization group approach, was set up 
in the invited lectures of Perry [Q], where the leading order calculation (deriving TYq) 
was completed, and the two calculations of this dissertation were mentioned as future 
prospective calculations — therein lies the historical motivation. 

Why is a calculation of the Lamb shift in hydrogen — which at the level of detail 
found in this dissertation was largely completed by Bethe in 1947 — or the ground 
state spin splitting in positronium — which through order a'^ was calculated by Ferrell 
over 40 years ago — of any real interest today? While completing such calculations 
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using new techniques may be very interesting for formal and academic reasons, our 
primary motivation is to lay groundwork for precision bound-state calculations in 
QCD. These calculations provide an excellent pedagogical tool for illustrating light- 
front hamiltonian techniques, which are not widely known; but more importantly, it 
presents three of the central dynamical and computational problems that we must 
face to make these techniques useful for solving QCD: How does a constituent picture 
emerge in a gauge field theory? How do bound-state energy scales emerge nonper- 
turbatively? How does rotational symmetry emerge in a nonperturbative light-front 
calculation? These questions can be answered directly in QED, as this dissertation 
shows. In QCD, the answers clearly change, but the overall computational framework 
does not, and thus an analytic understanding of the framework is essential. And, as 
already mentioned, QED allows this analytic understanding. 

An outhne of this dissertation follows. There are six chapters of which this is the 
first. The final chapter contains a general discussion and summary. In Chapter 2 we 
introduce light-front field theory: A pedagogical introduction to hght-front coordi- 
nates is given, and we present a simple tree-level example illustrating the vanishing of 
vacuum mixings in light-front perturbation theory; scalar theory is used to illustrate 
the division between kinematic and dynamic Poincare generators in light-front field 
theory; a derivation of a canonical QED Hamiltonian in the light-cone gauge is given. 
In Chapter 3 wc give an overview of the light-front hamiltonian bound-state problem 
and then discuss three renormalization group transformations in hamiltonian theory 
that are of interest; in the final section we discuss renormalization in light- front field 
theory which leads us to introduce "coupling coherence" which is elucidated with a 
simple one-loop example in coupled scalar 0{2) theory. Chapters 4 and 5 contain the 
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heart of this dissertation where the aforementioned calculations in positronium and 
hydrogen respectively are presented. In Chapter 4 we also present a simpler method 
of calculating the spin splitting, which may turn out to be useful in carrying out 
future higher-order calculations. 
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CHAPTER 2 



LIGHT-FRONT FIELD THEORY AND CANONICAL QED 

HAMILTONIAN 



In this chapter a pedagogical introduction to hght-front field theory, including 
a simple example to illustrate the vanishing of vacuum mixings in amplitudes and 
a discussion of the ten Poincare generators, is presented. Then a derivation of a 
canonical QED Hamiltonian is given. 

2.1 Light-front field theory 

Light-front coordinates (also called front form, null plane, or /zg'/it-cone coordinates — 
the usual coordinates are called instant form or equal-time) were first presented by 
Dirac in 1949 in his pursuit of alternative forms of relativistic dynamics that combine 
"the restricted principle of relativity with the Hamiltonian formulation of dynam- 
ics 0, 1^." Recent interest in light-front coordinates continues this pursuit in mainly 
two arenas: the low-energy nonperturbative bound-state problem of QCD, where 
light-front coordinates may allow a simpler vacuum structure (replacing the vacuum 
structure with the appropriate effective interactions through renormalization) than 
with equal-time coordinates, and high-energy scattering observables where light-front 
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coordinates are the natural coordinates of the system. For an extensive hst of hght- 
front references through the early 90s see 0. We must apologize for the inadequate 
recent references given to this diverse and active area of research, but for a fairly 
complete list see the following recent reviews [rU| and references within. 



An introduction to light-front field theory will now be given.Q An "ET" label 
implies an equal-time vector; a "LF" label implies a light-front vector. Light-front 
coordinates are defined by 

• ^LF ~ ^BT + ^BT (2-1) 

• ^LF — ^ST ~ ^FT (2-2) 

• "^LF ~ ■^ETy ^ = 1) 2 (2-3) 

• ^LF — (^LF) ^LF) ^Lf) (2-4) 

• ^FT ~ (^FT) ^FT) ^Ft)- (2-5) 

Carrying around these "LF" and "ET" labels is tedious, so following convention (an 
obvious convention given the signs on the right-hand-side of the top two equations 
above), we define 

x+ = xIp, (2.6) 
X = x^p , (2-7) 

and then drop the labels. There is no notational ambiguity b In 
this introduction we will keep the labels if possible confusion may arise otherwise. 
More formally the above is written 

^FF — "^^u ^ET • (2-8) 
^For some other works with nice introductions see for example [Tl| and references within. 
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Actually A^^^ is defined to be the same for all 4-vectors, but here we just show the 
transformation on the coordinate x^. Note that x^p and not related by a 

Lorentz transformation since det(y4) = —2, while a Lorentz transformation L must 
satisfy det(L) = ±1. The next step is the requirement that all Lorentz scalars are 
equivalent — done by adjusting the metric tensor. For the x^x^ scalar we have 

dj^J ^ET^ET = djlu ^LF^LF ■ (2-9) 

Since A = A^"^ (but note A ^ A^^) this simply implies 

9aP = ^Ya 9^1 ^Tp ^ 9lJwer = A ^ ■ gfo^^r ' ^ ^ (2-10) 

and 

9fF = A\ g^^p A^^ g^r = A ■ glT' " A , (2.11) 

where the terms on the far right are written in convenient matrix notation. Given the 
standard equal-time metric tensor, gp'rp = gj^J = (1,-1, —1, —1), we follow convention 
and drop the LF labels, take — > + and 3 — > — and end up withQ 

1 = ^ = ^ = 2(?+- = 2(7-+ = -g'' = -g^' = -gu = -g22 ■ (2.12) 

The components of the light-front metric tensor not mentioned are zero. Note that 
these factors of two in the metric tensor lead to factors of two in other places like 

1 11 

d^x = -dx^dx~(fx'^ , x_ = g-+x^ = -{x^ + x^) = -^^^o ~ ^3) > etc. 

Conventionally, x~^ is chosen to be the light-front time coordinate. Thus is the 
light-front longitudinal space coordinate. Also, from the p ■ x scalar, 

Pfj,x^ = g-^-^p~x~^ + g^j^p'^x~ — = -p'x'^ + -p'^x' — p*x* , 

^Note that this replacement '0 — > +' and '3 — > — ' appHes whether it be an upper or lower 
index — simple example: x^j^p = x+ and = x+ = = x^ /2. 
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we see that p~ is the hght-front energy coordinate, and is the hght-front longitu- 
dinal momentum coordinate. The light-front dispersion relation for an on mass-shell 
particle is interesting: 

= ^ p- = ^ . 2.13 

There is no square-root (compare p° = ±v^pM-~n7? in equal-time) and small longi- 
tudinal momentum p"*" implies high energy (except for a set of measure zero for a 
massless particle, i.e. m = p"*" = 0). 

All trajectories in the forward light-cone in light-front coordinates have p'^ > 0. 
The Lorentz-invariant 3-momentum integral — which is the method of summing over 
particle momenta whenever required — shows this (it also reminds us that particle 
lines in hamiltonian diagrams are on mass-shell): 

27r6{p'~my{p')f{p)= f -^f{p) 



cP'p^dp^9{p 



Especially note this last definition of , which is a shorthand used in the dissertation 
often. 



The consequences of p"*" > are illustrated by a simple example [T^ . Consider the 
Feynman diagram of Figure |2.1| for the tree-level annihilation amplitude in cfr' theory. 



This amplitude is a Lorentz scalar (and thus light-front and equal-time field theory 
should give the same result for this amplitude) given by 

^ (2.14) 



s — m? 



where 



S=[P1+P2Y{PI+P2)^- (2.15) 
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P1 P3 




Figure 2.1: Feynman diagram for the tree-level annihilation amplitude in 0^ theory. 




Figure 2.2: (a) Time-ordered diagram for the tree-level annihilation amphtude in 
(f)^ theory, (b) Time-ordered diagram containing vacuum mixings. It vanishes in 
light-front field theory. 
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In a hamiltonian approach all time orderings must be included which leads to the 



two diagrams of Figure p.2| . In equal-time field theory these diagrams contribute 

2E(pi + P2) [^(pi) + E{p2) - ^(pi + P2)] ^^'^^^ 

and 

(2.17) 



2E(pi + P2) [-^(pi) - E{p2) - E{pi + P2)] 
respectively, where -E'(p) = ^/p^~+~?7?. Summing Eqs. (|2.16|) and (|2.17| ) gives ^ 



In light-front field theory Eq. ( |2.17|) vanishes and Eq. (|2.16| ) becomes 



^-^TCT = ' (2-18) 



{pi+P2Y[P-{pi) + p-{P2)-P-{pi+P2)\ 'P+ \r±l+± - ^±1+Ilf.] S-im? 



•p+ p+ 

where P~{p) = ^ p+ , the 3-momentum V = pi + p2, and s is the same variable 
as defined in Eq. (|2.15| ). Thus the positivity of p~^ makes the diagram with vacuum 



mixings — as in Figure |2.2| b — vanish, and all the amplitude is in the one time-ordering 
alone. This observation alone sparked much of the initial interest in light-front field 
theory 0, and continues to be a topic of focus today. Today this topic deals with a 
subtlety that will not be discussed much in this dissertation but is intensely studied by 
the practitioners of light-front field theory; this subtlety deals with p"*" = [|T^ (the so 
called zero- modes) or p+ — > |T^ (renormalization group approaches), and whether 
or not there is any nonperturbative physics hidden in this sector of the theory. From 
the example just shown it is clear that at least perturbatively the zero-modes do not 
contribute to the amplitudes. We will drop zero-modes initially and replace their 
effects through renormalization counterterms fixed by coupling coherence. Coupling 



coherence is explained in Section p.6| , and for a discussion on "p"*" renormalization 
group" see Subsection p. 6.1 . 
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To complete this light-front introduction we will discuss the ten generators of the 
Poincare group and how the division between the kinematic and dynamic operators 
is different in equal-time and light-front field theory. The discussion is not intended 
to be complete. Transverse and longitudinal surface terms are just dropped, and a 
flavor of the subject is given. For further reading and references consult |T3 . 



The Poincare generators are constructed from the stress tensor T^*^ — for a concrete 
example take scalar theory 

dC 

T"" = —— - g'^'C = df^cpd^cP - g^'^'C . (2.19) 

On an equal-time surface the generators are expressed as 

p/^ = y (fxT^^" , J^"" = J d^xix^T^" - x^T°^) , (2.20) 

with the usual definition of boosts and rotations as = and J*-' = eijkJ^ 
respectively.^ Note that (P*, J*) are kinematic and (P°, i^*) are dynamic generators 
as can be seen by the fact that the kinematic ones do not contain the hamiltonian 
density and therefore do not shift states off the initial surface = 0. 
On a light-front surface the generators are expressed as 

P^ = ^J d^x^dx-T^^" , J^"" = \j d^x^dx-ix^'T^^ - x^T+>') , (2.21) 

where the boosts are = —^J^^ and = J+*, and the rotations are = J^^ 
and = J^* — see below for the relation between the equal-time and light-front 
generators. Note that (P+, P\E\J^) are kinematic as can seen by the fact that they 

do not contain the hamiltonian density and therefore do not shift states off 

^£123 = 1- Roman indices i, j, ... take on the values 1, 2 and 3 in equal-time discussions but only 
the values 1 and 2 for light-front discussions. 
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the initial surface = 0. The transverse rotation generators are dynamic as can 
be simply seen by the fact that they contain the hamiltonian density and thus shift 
states off the initial surface. The longitudinal boost generator = —\J^~ appears 
to be dynamic since it contains the hamiltonian density; to see that it is kinematic it 
is important to note that the initial surface is x"*" = (if it were x"*" = constant, then 
would move states off the initial surface). A standard Lorentz transformation 
shows that the time coordinates in two frames boosted by are related by 

= x+exp(-cj) , (2.22) 

where uj is the rapidity of the transformation [recall v = tanh(u;) where v is the 
relative speed of the two frames] . So a;"*" = in the one frame is seen as x^' = in the 
other frame: is kinematic.^ This point will be discussed with interactions below. 
Note that should perhaps be called a longitudinal scaling generator instead of a 
boost generator since the action of in momentum space on all particle momenta 
is 

p+ ^ p+ exp(-^) , ^ , (2.23) 
and the exact Hamiltonian of the theory is transformed as 

P — > P exp(cj) . (2.24) 

The relation between the equal-time and light-front generators is easily found by 
using A^j^ of Eq. ( ^.81 ). is a four- vector, so we have simply 

^LF = ^^v^ET 5 (2.25) 

^This simple example shows that even in free field theory in order for to be kinematic, x'^ = 
must be the choice for the initial surface. 
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which gives 



p+ — pO pO I p3 p— — p3 pO p3 

r — JTji^p — r^rp -t- r^np , r- _ rj^p — r^j, r^rp , 

= -Plf ~ -^ET 5 = -^LF ~ ^ET ■ (2.26) 



For the boosts and rotations, which form a tensor, we have 



T/J-u _ AfM jal3 Au ^ jupper _ a jupper a icy 

"^LF — a "^ET ^ f3 ^ LF — ^ ' ET ' ^ J l^-^ ' ) 

where the relation on the right is written in convenient matrix notation. If the lower 
indices are desired we can just use the light-front metric tensor (it keeps track of all 
the factors of 2) which in convenient matrix notation gives 

jLF _ LF jupper LF /q qoN 

'^lower ~ yiower ' LF ' iJlower ■ \^.^0) 

Recalhng footnote ^ and our conventions below Eqs. (|2.20|) and ( p. 211) , this gives 



^ — "^LF ~ •^ET ~^ ^ET ' ^ — "^LF ~ '-^ FT ~^ ^ FT ) _ J^p — 2Kpj, , 

= J^p = —JpT + Kpp , = jf^p = Jpp + K%rp , J'^ = jY'p = Jprp (2.29) 



for the upper light-front indices, and 

jLF _ '^ET ^ET jLF _ ET ^ FT jLF _ ^ FT 
■^01 ~ 2 ' "'02 ~ 2 ' "^03 ~ 2 ' 

jLF _ "^ET ^ET jLF _ FT ^ FT jLF _ j3 (c) nrW 

^31 ~ o ' 32 ~ o ' 12 ~ "^ET I^Z.OUJ 



for the lower light-front indices. To insure that all the factors of two are right, a nice 
check is 

Tr ( J-^P- ■ J,,^e.) = Jip J^f = J^T = 2 (K^ - J2) . (2.31) 
This holds true for the above relations. 
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To close this section we work out the scalar example more explicitly as an illus- 
trative example of how the interactions enter the generators. For fermion and vector 
examples see [0. Recall the form of and J^^ from Eq. ( p.21| ), where the stress 



tensor T^^^ was written in Eq. ( p^.l9| ). For concreteness, say the lagrangian density is 

1 Tn?' 
C = -d,<P d^(t> - ^ - ^inM) ■ (2.32) 

Thus we have 

T++ = (9+0)2 , (2.33) 
T+* = 9+0 (9> (2.34) 

and 

T+- = d+(j) d'(f)-2C 

= (9,0)2 + mV + 27^,,i(0) . (2.35) 

For scalar theory the kinematic generators are (note there are seven which is one 
more than in equal-time field theory) 

P+ = ^1 dVdx- (9+0)2 , (2.36) 
1 r .2 , , 



P" = -j d'x^dx-d+(f) 9> , (2.37) 
1 r 

= -J (fx^dx- (x^9+0 9^0 - ^29+0 9V) , (2.38) 

W = x+P' - - / d^x^dx-x\d+^y , (2.39) 



2 

X+P- 



+ d^x^dx-x~{d+(f)f . (2.40) 



This last equation makes it clear why is kinematic only for fields initialized at 
a:+ = 0; is the exact Hamiltonian of the theory and its job is to move states off 
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the initial surface, but at x"*" = this term in vanishes.F| The dynamic generators 



are 



= \j d^x^dx- [{di(pf + m202] + j (fx^dx^ HinM) , (2.41) 
F' = ^J (fx^dx- {x-d+<p d'(j) - x' [{di(f)f + m^^^j j 

- J d^x^dx'x' TiinM) ■ (2-42) 

Note that the stress tensor is symmetric and conserved: T^^^ = T^^ and d^T^^ = 
0.0 Thus the ten Poincare generators are independent of time x'^. For example, we 
will verify this for K^. Explicitly taking a derivative of Eq. ( ^.4(J| ) gives 
dK^ P- x+dP- If,, , . c/T+^ 



, , +- d^x^dx'x-——. 2.43 
dx+ 2 2 dx+ AJ dx+ ^ ^ 

The second term on the right is zero by 

1 r dP~ 

= - / d^x^dx~daT^~ = —— + surface terms , (2.44) 

2 J dx+ 

and the third term can be rewritten using 

= 1 y d^x^dx-x-d^J^'^ = ^J d^x^dx~x-{d+T++ + d^T'^) 

+ transverse surface terms ; (2.45) 

performing an integration by parts over x~, this becomes 
1 r dT^'^ 1 f 

- I d^x'^dx~x~ — — — = - / d'^x'^dx'T^^ + surface terms . (2.46) 
4 J dx+ A J 



Inserting Eqs. ( p.44| ) and ( |2.46| ) into Eq. ( |2.43| ) gives 
dK^ 1 r 

= -^j d^x^dx- (t+- - r-+) + surface terms = , (2.47) 

^Interestingly, note that is independent of time for all time as detailed below, but only if the 
initial surface is dialed precisely to = 0, does it become kinematic. 

^ The stress tensor is conserved only if the equations of motion a(d'4>) ~ satisfied, as is 

easily verified. 

14 



where in this last step surface termsQ are dropped and the fact that the stress tensor 
is symmetric is used. Similar algebra can be used to show that all ten Poincare 
generators are independent of time x~^. As a final remark note that the — ^"^^ term 
m was essential to show that is independent of time for all times. 

To review, P and are the "Hamiltonians" of the system, and K^, E\ J^, 
and are the kinematic generators of Lorentz transformations of the system. 
The fact that boosts are kinematic leads to a simple change of coordinates^ making 
this boost invariance manifest and the total momentum of the system is seen to 
factor completely out of the Schrodinger equation, lowering the dimensionality of the 
equation and making its analysis simpler (analogous to the nonrelativistic hydrogen 
problem in equal-time coordinates where the center-of-mass and relative coordinates 
factor in the Schrodinger equation). Note that the eigenvalue of (with = 0) is 
the helicity of the respective state — in this dissertation, the helicity of the electron 
will be written as s/2, where s = ±1. The fact that the rotation generators and P^ 
are dynamic makes classifying their states complicated. However, they are symmetry 
operators: [P^, P*] = 0, thus in principle they should lead to simplification in the 
diagonalization of P^. In practice, no such simplification has been found to date, 
thus we do not discuss them further in this dissertation. 

2.2 Canonical QED Hamiltonian 

A derivation of the canonical Hamiltonian of hydrogen treating the proton as a 

point particle is presented. For positronium neglect all terms that contain a proton 
^The point here is not to discuss surface terms, but to note the already present non-trivial 



cancelation in Eq. (2.47) 



^These coordinates are called Jacobi coordinates. More on this later. 
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field. A summary of tlie derivation and results will be given, and then the details of 
the derivation will follow. 

Starting with the QED lagrangian density for the electron, proton, and photon 
system (e > 0) 

= -^F^uF'^'^ + i^,it ^ + e^-m)ilJe + i^pii ^-e4-mp)i;p , (2.48) 

in a fixed gauge, A'^ = 0,0 the constrained degrees of freedom are removed explicitly, 
producing a canonical Hamiltonian H. m and rrip are the renormalized masses of 
the electron and proton respectively. ^ = ^4^7^. Details of the derivation follow 
below. For our 7-matrices we use the two-component representation chosen by Zhang 



and Harindranath [0. The field operator expansions and light-front conventions are 
summarized in Appendix |^. The resulting canonical Hamiltonian is divided into a 
free and interacting part 

H = h + v. (2.49) 

h is the free part given by 
h= fT. UipMp) + B\{v)BM ^— ^ + ^\{p)o^s{p) ^) ,(2.50) 

hs\ P^ P^ P^ j 

plus the anti-fermions. The notation for our free spectrum is h\i) = ei\i) with 
J2i = I5 where the sum over i implies a sum over all Fock sectors, momenta, and 
spin. Recall, we use the shorthand Jp = J ^iQ^apl^^^ ■ v is the interacting part given 
by 

V = j d^x^dx-Hint , (2.51) 

^This derivation will not include a discussion of the gauge field zero-modes. Initially we drop 
zero-modes. Their effects are conjectured to be replaceable by effective interactions that are fixed by 
coupling coherence. Coupling coherence is explained in Chapter^. For a treatment that incorporates 
these gauge field zero-modes from the start in QED see and references within. 
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where 

"^int 'Hee-y ~l~ '^ppy ~l~ "^6677 ~l~ '^ppy-y ~l~ "^7— inst ; (2.52) 



and 




Ke7 = eel <i -2(^ ■ A^) + a ■ + ^- a ■ } , (2.53) 

H.7 ^ 2(fi . A-) - . . A-^^^ - 



.9^ 



Ke77 = jela ■ A^^ia ■ A^Q^ , (2.55) 

^PP77 = -^e^ jeja ■ A^^ia ■ A%)^ , (2.56) 

^7-inst = '^^777772'^^ ■ (2.57) 

2 (9+) 

J"*" = 2e {^l^p — ^l^e) and a* are the standard SU{2) Pauh matrices, i = 1,2 only; 



e.g., (T ■ 9^ = cr*9* = (^^(^ gfr) + o"^(— gfj). The dynamical fields are A\ and 
^p, the transverse photon and two-component electron and proton fields respectively. 
For the relation between ijj and and a summary of our light-front conventions see 
Appendix 0. 

In this dissertation the formal expression for the canonical Hamiltonian given by 
Eqs. ( |2.49D -( p.57D is of little practical use. Of more practical use are matrix elements 
of the canonical Hamiltonian H in the free basis of h. The matrix elements used in 



this dissertation, as derived by Allen ]19[, are given in Appendix pj. 

Now the details of the derivation of the canonical Hamiltonian are presented. 
Given £ of Eq. (|2.48|) , the equations of motion are 

df.F^"' = r, (2.58) 

(i ^ + e 4 - m)i)e = , (2.59) 
17 



{i ^ — e Jji — nipjipp — , 



(2.60) 



where J'^ — e (^jjpj'^ijjp — i/j^'j^i/je^ ■ The physical gauge A'^ = is chosen and the 
projection operators A_|_ and A_ are inserted into the equations of motion. Note 
■0- = A-'^ and ■0+ = A^'^. Three of the equations are seen to be constraint equations: 



^d+d+A- + d'd+A' = J+ , 



e+ 



and 



(2.61) 
(2.62) 

(2.63) 



The fact that these are constraints can be seen from the fact that no time derivatives 
d~ appear. Note a* = 7^7*. Inverting the space derivative 9+ gives 



—2 (9* 



and 



(2.64) 
(2.65) 

(2.66) 



To get a flavor of a position space representation of the inverse longitudinal deriva- 
tive, note that it could in principle be defined as follows (see below for what we 
actually use), where j{x) is an arbitrary field: 



/ 1 \ ^ 1 



^ -y )f{y ) + 9i , 



-y \f{y )+92 + x gs, 



_d_ 
dx' 
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d-eix -y ) = 25{x -y ) , 

e{x) = e{x) - e{-x) . 

X"*" and X"*" are implicitly in the arguments of the f{x)s. gi, g2 and gs are arbitrary 
fields independent of x~ . For a discussion on these boundary terms see [^. Notice 
that this inverse longitudinal derivative is non-local. 

In practice, we define the inverse longitudinal derivative in momentum space. We 
explicitly put the momentum representation of the field operators into the respec- 
tive terms of the Hamiltonian, multiply the fields out explicitly, and then replace 
the inverse derivatives by appropriate factors of longitudinal momentum with the re- 
striction \p'^\/V^ > e = 0+ [P^ is the total longitudinal momentum of the physical 
state of interest]. The absolute value sign on \p~^\ is required for the instantaneous 
interactions. For example, a product of two fields gives 

exp[— ^(p — k) ■ x] — > — —6 (\p^ — k^\ — eV^^ exp[—i{p — k) ■ x] . (2.67) 



id~^ p+ — A;+ 

As a final note, we work in continuum field theory and drop all surface terms initially. 
If these terms are required for the Hamiltonian to run coherently, they are conjectured 
to arise through the process of renormalization. 

The dynamical degrees of freedom are A* , '?/'e+ and ipp^ . The canonical hamiltonian 
density is defined in terms of these dynamical degrees of freedom 

Taking these derivatives of the lagrangian density and combining terms, the hamil- 
tonian density takes the following simple form 

n = + ^pUd^^e- + ^Ud^i^v- - 2 (^^) + ' (2-69) 
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where the constraints of Eqs. ( |2.64|) - (|2.66|) are assumed to be satisfied. In our 7- 
matrix representation, only two of the components of the 4-spinors il)e+ and ipp^ are 
nonzero. Writing these as the 2-spinors and respectively,^^ and inserting the 
constraints of Eqs. ( |2.64| )-( |^.66| ), H takes on the form written earlier in Eqs. ( p^.49| )- 
(^.57]), where surface terms such as in 



j dPx^dx- (^^^) (^^^) = ~/ d^x^dx-J+ (^-^^ J++surface terms , (2.70) 
are dropped. 



^'^ In other words, we are defining tpi^^ — A+'f/^e = ('i/'ei) "062, 0, 0) = ('fe,0), with an analogous 
definition for the dynamical proton field. 
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CHAPTER 3 



LIGHT-FRONT HAMILTONIAN BOUND-STATE 

PROBLEM 



In this chapter, an outhne of the general computational strategy we employ is 
given. This includes a derivation of the effective Hamiltonian and a discussion of our 
diagonalization procedure. We discuss three hamiltonian transformations which are of 
interest ,Q give a general discussion of our renormalization program, and close with two 
simple examples: one introducing the standard renormalization group terminology 
through transverse and longitudinal renormalization group discussions, and the other 
elucidating coupling coherence. 

3.1 General computational strategy 

First, we give a brief overview. See the later sections of this chapter for the 

appropriate references and/or if the terminology is new to the reader. We use the 

renormalization group to produce a regulated effective Hamiltonian H\, where A is 

the final cutoff and renormalization is required to remove cutoff dependence from 

all physical quantities. This renormalization group transformation is run from the 

bare (initial) cutoff A down to the effective (final) cutoff A, with A ^ A. At this 
^^A fourth transformation is mentioned in a footnote. 
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point we have a regulated effective Hamiltonian that contains a finite number of 
independent^ relevant and marginal couplings — in the transverse renormalization 
group sense which is implied throughout this dissertation except for a brief example 
in Subsection |3.6.1| — and an infinite number of independent irrelevant coupling^ as 
would occur in any cutoff theory; also, there are an infinite number of new dependent 
relevant, marginal and irrelevant couplings. Note that we are speaking in terms of 
dimensionless couplings (with the appropriate power of the cutoff explicitly factored 
out). 

The reader may be concerned with the fact that there are an infinite number 
of dependent relevant and marginal couplings in the final Hamiltonian. This occurs 
because (i) longitudinal locality is absent and (ii) rotational symmetry is not manifest 
in light-front field theory. Predictability is not lost however because we require these 
dependent couplings to run coherently with the independent relevant and marginal 
couplings, that is we do not allow them to have explicit cntoS dependence. These extra 
conditions placed on the theory are called 'coupling coherence'. Coupling coherence 
fixes these extra counterterms, and the symmetries broken by the regulator (most 
noteworthy for us: gauge, chiral and rotational invariance) are restored in the final 

physical solution. These conjectures have been found to be true in all examples 

^^A coupling required in the theory when the symmetries of interest are maintained by the reg- 
ulator will be called 'independent.' Additional couplings required to restore symmetries broken by 
the regulator will be called 'dependent.' The reasoning behind this terminology is made clear below. 
Note that by 'symmetries' we include 'hidden symmetries': fixed relations between couplings in the 
broken phase due to an underlying symmetry of the theory which is however not maintained by the 
vacuum. 

Since A <C A, the independent irrelevant couplings become functions of only the independent 
relevant and marginal couplings, and are exponentially insensitive to their boundary conditions set 
at the bare cutoff A (universality). 
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worked to date; see Section |3.6.2| for a simple one-loop example and for references to 
further examples. 

This complicated effective Hamiltonian cannot be directly diagonalized, and since 
we want to solve bound-state problems we cannot solve it using perturbation theory. 
The second step is to approximate the full effective Hamiltonian using 

Hx = no + {Hx -Ho) = Ho + V , (3.1) 

where Tio is an approximation that can be solved nonperturbatively — for QED analyt- 
ically, which is one of the primary motivations for studying QED — and V is treated in 
bound-state perturbation theory (BSPT). The test of Tig is whether BSPT converges 
or not. 

We now formulate the following questions. Is there a range of scales A for which 
Tio does not require particle emission and absorption? What are the few-body inter- 
actions in 7^0 that generate the correct nonperturbative bound-state energy scales? 
Is there a few-body realization of rotational invariance; and if not, how does rota- 
tional symmetry emerge in BSPT? We should emphasize that for our purposes we 
are primarily interested in answering these questions for low-lying bound-states, and 
refinements may be essential to discuss highly excited states or bound-state scattering. 

It is essential that A,0 which governs the degree to which states are resolved, be 
adjusted to obtain a constituent approximation. If A is kept large with respect to 

all mass scales in the problem, arbitrarily large numbers of constituents are required 

^''The tilde on A is a useful notation since we want to discuss a regulator that restricts binding 
energy scales of a typical atom. Since we use different transformations for the positronium and 
hydrogen calculations the relation between A and A changes accordingly: Xhydrogen = ^ 2{rT+m^'')^ 
where a Bloch transformation is used and XposUronium = 2{m+m) ^^^^^ ^ G-W transformation is 
used. The difference (besides the trivial nip < — > m difference) arises simply because G-W regulates 
energy differences whereas Bloch regulates energies. 
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in the states because constituent substructure is resolved. A constituent picture can 
emerge if high free-energy states couple perturbatively to the low free-energy states 
that dominate the low-lying bound-states. In this case the cutoff can be lowered until 
it approaches the nonperturbative bound-state energy scale and perturbative renor- 
malization may be employed to approximate the effective Hamiltonian. In QED a 
simplification occurs because the system is nonrelativistic. For example in positron- 



lum 

= \ + 0{a), (3.2) 

^2^2*1^^ = 0{a) , (3.3) 
m 

a < 1 , (3.4) 

where m is the renormalized electron mass and is the total longitudinal momentum 
of positronium. Also, the dominant photon wavelength that couples to the bound- 
state is of order the size of the bound-state: wavelength ~ l/{ma). There is a natural 
gap in the system between the valence and valence plus one photon states and this 
above-mentioned range into which the cutoff must be lowered is 

ma^ <^ A <C ma , (3.5) 
where m is the electron mass and the relation between A and A is explained in foot- 



note TA. If the cutoff is lowered to this range, hydrogen and positronium bound-states 
are well approximated by electron-proton and electron-positron states respectively, 
and photons and pairs can be included perturbatively. 

It is an oversimplification to say the constituent picture emerges because the 
QED coupling constant is very small. Photons are massless, and regardless of how 
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small a is, one must in principle use nearly degenerate bound-state perturbation 
theory that includes extremely low energy photons nonperturbatively. This is not 
required in practice, because the Coulomb interaction which sets the important energy 
scales for the problem produces neutral bound-states from which long wavelength 
photons effectively decouple. Because of this, even though arbitrarily small energy 
denominators are encountered in BSPT due to mixing of valence bound-states and 
states including extra photons, BSPT can converge because emission and absorption 
matrix elements vanish sufficiently rapidly. A nice example of this is seen in the Lamb 
shift calculation of Chapter |^. 

The well-known answer to the second question above is the two-body Coulomb 
interaction sets the nonperturbative energy and momentum scales appropriate for 
QED. We have already used the results of the Bohr scaling analysis that reveals the 
bound-state momenta scale as p ~ ma and the energy scales s& E mo? . As a result 
the dominant photon momenta are also of order ma, and the corresponding photon 
energies are of order ma. This is what makes it possible to use renormalization to 
replace photons with effective interactions. The dominant photon energy scale is much 
greater than the bound-state energy scale, so that A can be perturbatively lowered 
into the window in Eq. ( |3.5|) and photons are not required in the state explicitly, but 
only implicitly through effective interactions. 

Finally we discuss rotational invariance in a light-front approach. In light-front 
field theory, boost invariance is kinematic, but rotations about transverse axes in- 
volve interactions. Thus rotational invariance is not manifest and all cutoffs violate 
rotational invariance in light-front field theories. In QED it is easy to see how coun- 
terterms in }i\ arise during renormalization that repair this symmetry perturbatively; 
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however, the issue of nonperturbative rotational symmetry is potentially much more 
complicated. We first discuss leading order BSPT and then turn to higher orders. 

To leading order in a constituent picture we require a few-body realization of 
rotational symmetry. This is simple in non-relativistic systems, because Galilean 
rotations and boosts are both kinematic. In QED the constitutuent momenta in all 
low-lying bound-states are small, so a non-relativistic reduction can be used to derive 
TYo- Therefore to leading order in QED we can employ a non-relativistic realization 
of rotational invariance. 

At higher orders in BSPT rotational invariance will not be maintained unless cor- 
rections are regrouped. The guiding principle in this and all higher order calculations 
is to expand not in powers of V, but in powers of a and a ln(l/a). TIq should pro- 
vide the leading term in this expansion for BSPT to be well-behaved, and subsequent 
terms should emerge from finite orders of BSPT after appropriate regrouping. Powers 
of a appear through explicit dependence of interactions on a, and through the de- 
pendence of leading order eigenvalues and eigenstates on a introduced by interactions 
in 7^0- This second source of dependence can be estimated using the fact that mo- 
menta scale as ma in the bound-state wave functions. Of more interest for the Lamb 
shift calculation is the appearance of a ln(l/a), which is signaled by a divergence 
in unregulated bound-state perturbation theory. As has long been appreciated, such 
logarithms appear when the number of scales contributing to a correction diverges. 

Now we change gears a bit. There is an open question as to which transformation 
should be used to derive the renormalized Hamiltonian. We discuss three different 
transformations: G-W, Wegner, and Bloch.|^ The G-W transformation is unitary and 

""^^A fourth transformation of promise is currently under study by Harada, Okazaki and co-workers 
PH . They combine the simplicity of Feynman perturbation theory with a similarity transformation 
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was developed by Glazek and Wilson |22| to deal with the small energy denominator 
problem. This transformation will be used for the positronium calculation. An early 
application which used the G-W transformation WM was a weak-coupling treatment 



of QCD. The Wegner transformation P3 is also unitary and allows no small energy 



denominators — the elegance in the definition of the formalism is one of the appealing 
features of the Wegner transformation. The Wegner transformation was developed 
with condensed matter applications in mind, however it has recently been applied 
to QED |2^. Another transformation — used for the Lamb shift calculation — is the 



Bloch transformation ||2^. Recent applications in scalar field theory p6[ and QED 



27| , utilize the Bloch transformation. It is not a unitary transformation and may have 
small energy denominator problems, but it is the simplest transformation to apply in 
practice and the natural gap in QED between the valence and valence plus photons 
sectors avoids problems here. A final note is that the small energy denominator 
problem should not arise in general when solving for the low-lying spectrum which 
we discuss next. 

The small energy denominator problem is explained as follows. An effective Hamil- 
tonian with energy cutoff A has effective interactions in the low-energy space L at 
second order of the following generic form 

{L H^jf ^) ~ 2^ 7^ r- , (3.6) 

where \L) and \H) are in the low and high energy spaces respectively and Eh is above 
A (and below A) while El is below A. This effective interaction arises to replace 

the effects of couplings removed by lowering the cutoff from A to A. The small 

that moves the effects of couphngs between few and many-body states into effective interactions 
resulting in an effective Hamiltonian acting in the few-body space alone. 
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energy denominator problem arises because in general can go all the way up to 
the boundary A and Eh can go all the way down to the boundary A. This can lead 
to a vanishing denominator which is problematic unless the matrix elements vanish 
fast enough and lead to either a vanishing or finite result. Note that if El (which 
is not integrated over in J2h but rather is fixed by the wave function) is much less 
than the boundary A then there is no problem. Thus we want to keep A relatively 
high in comparison to the nonperturbative bound-state energy scale. But also as 
detailed above we want to keep A relatively low to suppress many-body mixings such 
as photon emission. In QED there is a range allowed for A that is high enough, so 
the small energy denominator problem is not encountered, and at the same time is 
not too high, and thus photon emission can be removed perturbatively. This range is 
written above in Eq. ( |3.5| ). In summary, even with the Bloch transformation, the small 
energy problem should not arise when solving for the low-lying spectrum because A 
can always be adjusted to be well above this level and \El — Eh\ ~ A or bigger, and 
can not vanish. In addition, in QED with the natural gap |-Eee7 — Eee\/ E^e ~ 137, 
we can still keep A low enough and perturbatively replace emission effects with few- 
body interactions in the valence sector alone. The small energy denominator problem 
should only begin to arise when for example calculating a scattering amplitude of 
energy near the boundary A. 

3.2 Step one: Derivation of Hx with G-W transformation 

A self-contained discussion on the derivation of the effective Hamiltonian Hx with 
the G-W transformation will now be given. The discussion will be general, and in 
particular will hold for both QED and QCD. 
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The derivation starts with the definition of a bare Hamiltonian H\ 



Ha = h + v^, (3.7) 

Va = /a^a . (3-8) 

^A = v + ^^A^ (3-9) 

H = h + v, (3.10) 

where h is the free Hamiltonian, H is the canonical Hamiltonian,Q is a regulat- 
ing function, and 6v^ are counterterms defined through the process of renormaliza- 
tion. The canonical Hamiltonian H is written in terms of renormalized parameters 
which implicitly depend on the renormalization scale /i as explained below. The 
counterterms are fixed by coupling coherence. In Section ^]6| coupling coherence 
will be explained further and also a discussion of the scale dependence of the theory 
will be given. 

The free Hamiltonian h is given by (using positronium as an illustration) 

h = /E r^t'"' [bliP)bs{p) + dlip)dsip)]+^al{p)a,{p)\ , (3.11) 

•J V o \ jJ U I 



h\€)=ei\{) , ^|z)(«| = l, (3.12) 



where the sum over i implies a sum over all Fock sectors and spins, and integrations 
over all momenta in the respective free states, m is the renormalized fermion mass. 
The regulating function is defined to act in the following way 



(^I/a^aIj) = fMM^A\j) = Imj^AH ^ (3-13) 

^^In the initial setup of the G-W transformation, A and A will be used as a shorthand for ^ 
and ^ respectively (in other words they are to be thought of as having the dimension of energy), 
where 7^+ is the total longitudinal momentum of the physical state of interest, and and have 
dimension (mass)^. 



^^See Section 2.2 for a canonical QED Hamiltonian. See for a canonical QCD Hamiltonian. 
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/,,^. = ^^(A-|A,,|) , A,, =£,-5,. (3.14) 

Note that this choice of a step function is not necessary and can lead to pathologies, 
however it is useful for doing analytic calculations. 

Next, a similarity transformation is defined that acts on H\ and restricts the 
energy widthsQin the effective Hamiltonian Hx to be below the final cutoff, A (A < A). 
This transformation allows recursion relationships to be set up for H\, which can be 
written in the following general form 

Hx = h + vx, (3.15) 
v\ = f\v\ , (3.16) 
v\ = Vx +Vx ^ , (3.17) 

where the superscripts imply the respective order in the canonical interaction v, which 
recall is written in terms of renormalized parameters. 

Now, starting with the above bare Hamiltonian, we will describe this procedure 
more explicitly. The similarity transformation is defined to act on a bare cutoff 
continuum Hamiltonian, H/^, in the following way: 

Hx = S{X,A)HaS\X,A) , (3.18) 
S{X,A)S\X,A) = S\X,A)S{X,A) = 1 . (3.19) 

This transformation is unitary, so and Hx have the same spectrum: 

Ha\^a) = E\^a), (3.20) 
S{X,A)H^S^{X,A)S{X,A)\^^) = ES{X,A)\^a) . (3.21) 

""^^The magnitude of the energy difference between the free states in a matrix element of a Hamil- 
tonian is defined to be its energy width. 
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Therefore, E is independent of the final cutoff A if an exact transformation is made. E 
is also independent of the bare cutoff A after the Hamiltonian has been renormalized. 

To put the equations in a differential framework, note that Eqs. ( p.l8| ) and ( p.l9| ) 
are equivalent to the following equation (proven below) 



d\ 

with 



[Hx,T,] , (3.22) 



S{X,A) = Texp(^l\x'dX'^ , (3.23) 

where Ty is the anti-hermitian (T| = —Tx) generator of energy width transformations, 
and T orders operators from left to right in order of increasing running cutoff scale 
A'. Eq. ( p.22| ) is a first order differential equation, thus one boundary condition must 



be specified to obtain its solution. This boundary condition is the bare Hamiltonian: 
Hx\^_,^ = H\. H\ is determined by coupling coherence. 

Now we prove that Eqs. ( |3.22| ) and ( |3.23| ) are equivalent to Eqs. ( |3.18| ) and ( |3.19| ). 



First note that to prove unitarity, Eq. ( p.l9|) , we need to know 



S^{X,A) = T+exp j\ydX'^ , (3.24) 

where we used the fact that Ta' is anti-hermitian, and orders operators from left 
to right in order of decreasing running cutoff scale A'. Now for the proof: Take a 
derivative of Hx and use unitarity giving 



dH 



X ^ ^^(^iJS}.s^^X,A)S{X,A)HjyS^{X,A) 



dX dX 

"'^^Note that the bare Hamihonian is not just the canonical Hamiltonian written in terms of bare 
parameters. The bare Hamiltonian must be adjusted through the process of renormalization until 
the conditions of coupling coherence are satisfied. 
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+ 5(A,A)ifA5t(A,A)5(A,A) 



rfS't(A,A) 



dS{X,A) 



S\\,A)Hx + HxS{\,A) 



d\ 
dX 



(3.25) 



Taking a derivative of the unitarity condition gives 

dS^{X,A) dS{X,A) 



and thus we have 



5(A,A)- 



dH. 



dX 



dX 



S\X,A) 



dX 



Hx,-^^%^S\X,A) 



Also, exphcitly taking a derivative of Eq. (|3.23| ) gives 



dS{X,A) 
^A 



exp 



TydX' 



-TxS{X,A) 



(3.26) 



(3.27) 



(3.28) 



which after using unitarity combined with Eq. (|3.27| ) imphes 

dH^ 



dX 



(3.29) 



as was to be shown. So, we must specify Tx and the running of Hx with the cutoff is 
completely defined, and it runs unitarily. 

To define Tx note that it is enough to specify how vx and h change with the 
running cutoff scale A. This is seen by writing out Eq. ( |3.29| ) more explicitly using 



Eq. ( CT) : 



^ + ^(/At^A) = [h,Tx] + [vx,Tx] 



(3.30) 



We solve this perturbatively in v, and note that h depends only on the renormalization 
scale /i (more on this later) but not on the running cutoff scale A. Also, we demand 
that Tx and vx do not contain any small energy denominators. Thus we define 



dh 

dx - °' 



(3.31) 
(3.32) 
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Eq. (|3.32|) is a choice such that Tx and consequently vx do not allow any small energy 
denominators. These additional constraints determine Tx and vx, which are given by 
the following equations 



[h, Tx] 

Vx 



_ dfx -J r , 

- fx[vx,Tx] , 



V + 6v. 



[vx',Tx']d\' 



(3.33) 
(3.34) 



where fx + fx = I ■ Eqs. (|3j|) and (^M ) follow from Eqs. (gjg)-(|3j|) and the 
boundary condition Hx\^^^ = H\. Now we solve Eqs. ( p.33| ) and (|3.34|) for Tx and 
Vx- Given Eqs. (p.l5|) and ( p.l6|) , we need to determine vx, and Hx is then known. 
The perturbative solution to Eqs. (|3.33| ) and (p.34|) is 



Vx 



_(1) , _{2) 

Vy. +V. + 



Tx = Tx' + Tx' + 



(3.35) 
(3.36) 
(3.37) 



where the superscripts imply the respective order in the canonical interaction, f , and 
these quantities are given by 



KTx 



(1) 



_(2) 



h,Tx 



_(3) 



h,Tx 



dfx 



dX ' 



_(2) dfx 



-F r (1) rr^W 

Vx^-fx[vx.Tx 



_(3) dfx 



dX \Vy , Ty 

Ix (K'.^; 



+ 

(2) 

A 



(2) ^(1) 

^A' . Ty 



+ 5v 



(3) 



+ 



(2) ^(1) 

^A ,Tx 



(3.38) 
(3.39) 
(3.40) 
(3.41) 
(3.42) 
(3.43) 
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A general form of these effective interactions is 

vx = -EW.)/ dX'[v'^\TS'] + Sv'^ , (3.44) 



for i = 2,3, ■ ■ ■, with f = v. 



Recall the general form of the effective Hamiltonian given by Eqs. ( p.l5|) -( ^.17| ). 



The explicit perturbative solution through third order now follows, v^^^ = v, and 
now we write the explicit form of the second order effective interaction . From 



Eq. ( |3.40D we obtain 



(AA) (AA)\ 

_ „. / 3ikj , 9jki \ , ,J^) 

Xij 



V 



E^^^^^.R^ + T^ (3.45) 



where gitf^ ^ T d\' fy,,^ . (3.46) 



,w. dfx> 
A dX 

(2) 

Svf^ij will be determined by requiring the conditions of coupling coherence to be 
satisfied.^ These previous equations are valid for an arbitrary regulating function 
fx- In this work we will use fxij = 6{X — |Ajj|) [a convenient choice for doing analytic 
calculations]. This gives 

fi'ifcf ^ = ifA-ik - fxik) Qikj , (3.47) 

Qikj = ^(|A,fc| -|Afc,-|) . (3.48) 

(3) 

For completeness, the explicit form of the third order effective interaction Vx is 
written in Appendix |C[ 



^°For the example of positronium's self-energy see Subsection |4.1.l| . See Subsection 3.6.2 for a 
coupled scalar theory example. 
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3.3 Step two: Diagonalization of Hx 

The second step in our approach is to solve for the spectrum of Hx- The Schrodinger 
equation for eigenstates of H\ is 

Hx\^.A'P)) = E,\^,,{V)), (3.49) 
which written out in all its many-body complexity is 

Y.^\Hx\j)m,^A'P)) = E,(^\^,^,{r)) , (3.50) 

3 

where the sum over j implies a sum over all Fock sectors and spins, and integrations 
over all momenta in the respective free states | j') . W labels all the quantum numbers 
of the state, and is discrete for bound states and continuous for scattering states. 

= — ^ , where for positronium for example the binding energy —Bn is defined 
by = (2m + B^^y. 'P' is the total momentum of the state of physical interest. 
Solving this eigenvalue equation exactly is not feasible, because all sectors are still 
coupled. For example, for positronium we have 

l*x,.(^)) = E = E \ee{if)){ee{if)\^,jV)) 

i i' 

+ E |ee7(^'))(ee7(^')|*.,.(P)) + E |eeee(z'))(eeee(z')|*,,,(P)) + • • • , (3.51) 

i' i' 

where the sum over i' implies a sum over all spins, and integrations over all momenta 
in the respective free states 

As already mentioned, we divide Hx into two pieces 

Hx = no + (Hx-no) = no + v , (3.52) 

diagonalize Hq exactly for these QED calculations, and calculate corrections to the 
spectrum of Hq in BSPT with V to some consistent prescribed order in a and 
Q;ln(l/Q;). 

35 



We close this section by writing the standard BSPT Raleigh- Schrodinger formulas. 
For simplicity, we write the formulas for the non-degenerate case p8[ : 



{no + V)\^,^,{V)) = E,\^iI,jV)) , (3.53) 
HolijAV)) = SJ^P.iV)) , (3.54) 
W> = I'PJV)) + E " + O(V') , (3.55) 



2 



(V'jv(7')|^jv(^)) 



C(V^) , (3.56) 



where P is the total three-momentum of the state and "N" labels all the quantum 
numbers of the respective state. These formulas will be used in Chapter ^ to solve 
for positronium's leading ground state spin splitting (actually degenerate BSPT is 
required here) , and in Chapter ^ to calculate the dominant part of the Lamb shift in 
hydrogen. Note that for the light-front case = — ^ and 8^^ = — ^ , where 
and are the exact and leading-order mass-squared respectively. 



3.4 Wegner transformation 

The Wegner transformation is unitary and thus follows the initial discussion on 
the G-W transformation in Section |3.2|. The Wegner transformation is defined by 



Tg = [Hg, and so the Hamiltonian evolves according to 

dHg 



ds 



Hs, [Hs, h]] , (3.57) 



where s is a cutoff of dimension 1/ (energy)^, which is obvious from the form of this 

defining equation. The free Hamiltonian h only depends on the renormalization scale 

^^Actually, Wegner uses Hd (the full diagonal part of the Hamiltonian at scale 's') instead of h 
(the free Hamiltonian). We choose the free Hamiltonian for its simplicity. 
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fx but not on the running cutoff s: |j = 0. See Section g]6| for further discussion on 
this point. Proceeding we define 

Vsij = exp{-sA^ij)vsij , (3.58) 

where h\i) = ei\i), Ajj = £j — (^l^sb) = Vgij and = h + Vs- After a little algebra 
we have 

—f^ = (^ik + ^ifc) "^sikVskj exp [-2sAikAjk] . (3.59) 
as ^ 

It was useful to note A^^- — Af^ — A|^ = —2AikAjk. Note that this is completely 
nonperturbative so far, however perturbatively it does start out at second order on 
the right-hand-side. This transformation has a simple nonperturbative form, but to 
date has only been solved perturbatively in field theory. We make the perturbative 
ansatz 

Vs = v's ^ + H , (3.60) 

where the superscript implies the order in the bare interaction Vg^ ■ At first order this 



gives 



—^=0, (3.61) 
as 

which implies 

%ii = ^s^ij , (3.62) 

where is the final cutoff. Note that > because of the dimension of the cutoffs, 
and the "no cutoff limit" is s„ — > 0. At second order we have 



B 

-(2) 



This exponential is easy to integrate with result 



Y{Aik + Ajk)Vs ikVs„kjQW[-'^sAikAjk] ■ (3.63) 



V 



(2) _ _ M , J 

~ 2 Y '''''''''''' K^^k ^iK 
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X [exp (-2s^ A,feA,fc) - exp (-2s^A,feA,fc)] . (3.64) 

We could continue this exercise to arbitrary order without too much apparent diffi- 
culty because exponentials are easy to integrate, but we will stop here since this trans- 
formation is not used in the specific QED examples of this dissertation. Basically, 
the Wegner transformation is a G-W transformation with the regulating function fxij 
chosen to be a Gaussian exp(— A^^/A^). 

3.5 Bloch transformation 

In this Section we discuss a derivation of an effective Hamiltonian via a Bloch 
transformation. As already mentioned, we use this transformation for our Lamb shift 
calculation. We use the Bloch transformation to separate the low and high energy 
scales of the problem and derive an effective Hamiltonian acting in the low-energy 
space alone with an identical low-energy spectrum to the bare Hamiltonian. We 



closely follow Section IV of |^ , where a discussion, including the original references, 
and derivation of a general effective Bloch Hamiltonian can be found. 
We start with a bare time-independent Schrodinger equation 

/^aI^a) = E|^a) . (3.65) 

Then projection operators onto the low- and high-energy spaces, Pl and Ph respec- 
tively, are defined. 



H 



v+ 
v+ 



h \ e\h 



v+ 



(3.66) 
(3.67) 
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Pl + Ph = ^ ^1 ' (3-68) 

where 6{x) is a step function. A and A are the bare and effective cutoffs respectively 
with A < A.0 V = (V'^, P"*") is the total momentum of the physical state of interest. 
h is the free Hamiltonian, which for hydrogen is written in Eq. ( p.50|) . 



An effective Hamiltonian acting in the low-energy space with an equivalent low- 
energy spectrum to H/<^ is sought. To proceed, a new operator 71 is defined that 
connects the Pl and Ph spaces 

PhI^a) =7^PL|^A) . (3.69) 



For a discussion and construction of TZ, see Eq. (4.4) and below in p9 |. 

This leads to the following time-independent Schrodinger equation for the effective 
Hamiltonian 

H^l^x) = ^|$a) , (3.70) 



where E is the same eigenvalue as in Eq. (|3.65|) , the state |$a) is a projection onto 



the low-energy space [with the same norm as |\E'a) of Eq. ( p.65|) ] 



I^a) = ^i + nmPL\^A) (3.71) 

and Hx is a hermitian effective Hamiltonian given by 

Hx = -j=1=={Pl + 7^t) Ha {Pl + TZ) , ^ . (3.72) 

Note that Hx acts in the low-energy space. In principle, all bare states |\E'a) that have 
support in the low-energy space have a corresponding eigenvalue given by solution to 



the eigenvalue equation (|3.7CI| ). 

22The t 
confusion 



^^The shorthands ^ ^4! — — > A and ^ — — > A are often used when it does not lead to 
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Defining H\ = h + v\, where h is the free field theoretic Hamiltonian and v\ are 
the bare interactions^ through fourth order in f a, the effective Hamihonian is given 
by 

{a\H,\b) = {a\h + VA\b) + lj:( + {a\vA\^){^Km 



Hi 



_^ 1 ^ / {a\vA\i){i\vA\j){j\vA\b) ^ {a\vA\i){i\vA\j){j\vA\b) \ 
_ 1 ^ / {a\vA\i){i\vA\c){c\vA\b) ^ {a\vA\c) {c\vA\i) {i\vA\b) \ 
+ {a\v[''\b), (3.73) 



where the fourth order terms are written in Appendix |C|. Aja = Si — Ea, with h\i) = 
ei\i). We are using \a), \b), ■■■to denote low energy states (states in Pl) and 
■ ■ ■ to denote high energy states (states in Ph). Note that the denominators 
are all Ea — £i — there are no Sa — £b or Si — Sj denominators from the same space. 
However, note that there are potentially problematic 

^mm (3.74) 

ZA^i . . . 



c.i 



and 

j^MMWM (3.75) 

type terms in the effective Hamiltonian beyond second order. See the already men- 
tioned Reference for a description of an arbitrary order (in perturbation theory) 



effective Bloch Hamiltonian and also for a convenient diagrammatic representation of 
the same. 

^^h is written in terms of renormalized parameters, and it is convenient to define va = v + Sva, 
where v is the canonical field theoretic interactions written in terms of renormalized parameters 
and S va ar e the c ounterterms that must be determined through the process of renormalization. See 
Eqs. (^^-(^) for the canonical Hamiltonian of the hydrogen system. 
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3.6 Renormalization issues in light-front field theory 

The renormalization concepts upon which this dissertation are based were intro- 



duced and elucidated by several examples in In |^ this work was continued 

with "A Renormalization Group Approach to Hamiltonian Light-Front Field The- 
ory." These studies were a continuation — generalizations to the light-front arena — of 



Wilson's seminal work on the modern renormalization group reviewed in |32| , p3 



with a nice simple introduction given in |3J]. Here we will give an overview of the 



concepts and elucidate them with a few examples. 

The effective Hamiltonian H\ is renormalized by requiring it to satisfy "coupling 
coherence." Simply put, our regulator breaks gauge and rotational invariance, and we 
need some principle in order to specify our counterterms in a unique and meaningful 
way, such that physical results are cutoff and renormalization scale independent, and 
manifest these symmetries. We demand that these additional counterterms required 
to restore the symmetries broken by the regulator run coherently with the "canon- 
ical renormalizable couplings" — that is we do not allow them to explicitly depend 
on the cutoff but only implicitly through the "canonical renormalizable couplings" 
dependences on the cutoff. This is coupling coherence.^ An example will be worked 
out in Subsection |3.6.2| . See the previously mentioned |^ for further examples, and 



note that "coupling coherence" was first formulated by Oehme, Zimmermann and 
co-workers in [^], where additional worked examples may be found.[^ 

The additional requirement of coupling coherence is necessary because without 

this assumption in a light-front effective Hamiltonian, there are an infinite number 

Actually, there is an additional requirement as explained below. 

Oehme, Zimmermann and co-workers use the term "reduction of couplings" instead of "coupling 
coherence." 

41 



of relevant and marginal couplings and predictability is lost.[^ This necessity was 
brought on by the fact that the only regulators we know how to employ in the non- 
perturbative bound-state problem in light-front field theory break Lorentz covariance 
and gauge invariance. Also, longitudinal locality (and perhaps transverse, but in this 
dissertation we assume transverse locality) is lost. The net result is that an infinite 
number of relevant, marginal and irrelevant couplings arise; that is, above and beyond 
the usual finite number of relevant and marginal couplings and the infinite number 
of irrelevant coupling^ as given by any cutoff effective field theory. These extra 
couplings are required to restore symmetries broken by the regulator or vacuum, and 
it is conjectured that coupling coherence restores these symmetries. For the discus- 



sion that follows, following Perry and Wilson |30|, call these extra couplings required 
to restore the symmetries of interest dependent couplings, and call the couplings re- 
quired even if the symmetries of interest are maintained by the regulator and vacuum 
independent couplings. Coupling coherence by construction fixes all the dependent 
couplings to be functions of only the independent relevant and marginal couplings 
near the lower cutoff A with A ^ A. 

In the notation of the G-W transformation, a coupling coherent Hamiltonian writ- 
ten in terms of dimensionless couplings for A ^ A satisfies 

Hx = S{\,A)HA{A,eA,mA,WA,c{eA,mA))S^{\,A) 

— > Ha{\, ex, mx, w{ex, mx), c{ex, nix)) , (3.76) 

with the additional requirement that all dependent couplings (represented by 'c' in 

the argument of the Hamiltonians) vanish when the independent marginal couplings 

^^The terminology "relevant," and "marginal," as well as "irrelevant," is explained below. 
^^Due to universality, predictability is not lost for this infinite set of couplings. 
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are set to zero. In Eq. ( p.76| ), ca and are independent dimensionless marginal 



and relevant couplings respectively, and wa represents the infinite set of independent 
dimensionless irrelevant couplings.^ c(eA, m\) represents the infinite set of dependent 
dimensionless relevant, marginal and irrelevant couplings. Coupling coherence fixes 
these couplings through 

dc{ex,mx) ^ dc{ex, mx) dex _^ dc{ex,mx) dirtx ,^ 
dX dex dX drux dX 

that is 

^ , (3.78) 

[note that this previous equation is implied in the notation of Eq. (|3.76 )] and as 
already mentioned, 

c{0,mx)^0. (3.79) 
The initial bare Hamiltonian if a does not satisfy Eq. (3.76), its form changes 



under the action of 5'(A,A). i^A rnust be adjusted until its form does not change. 
This "adjustment" is the process of renormalization. Coupling coherence is a highly 
non-trivial constraint on the theory and to date has only been solved perturbatively. 
In this dissertation, for QED, Eq. (|3.76| ) is solved to order e^, which turns out to be 



fairly simple because e does not run until order e . In Subsection |3.6.2| , solutions to 



Eq. ( |3.76D are obtained for the running of the marginal coupling through one loop in 



massless coupled 0(2) scalar theory. 

As far as the scale dependence goes, the canonical Hamiltonian H depends only 

on the renormalization scale /i which actually does not explicitly enter discussions 

^®As already mentioned, the independent irrelevant couplings flow to a function of the independent 
relevant and marginal couplings for A ^ A in the standard way according to the renormalization 
group. That is, the independent irrelevant couplings are automatically coherent and additional 
assumptions need not be made here. 
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of renormalization until including at least the running coupling e; the counterterms 
depend on A and this renormalization scale /i, but can not depend on the final cutoff 
A or the procedure is ill-defined. The effective Hamiltonian H\ depends only on A, 
however in order that it does not depend on the renormalization scale /i, e and m must 
implicitly depend on /i. If the renormalization procedure has been completed to some 
order in e, the limit A — > oo can be taken,0 and physical observables, in particular 
the spectrum, are independent of the remaining scales A and /i. The A-independence 
arises explicitly through the diagonalization process, while the /i-independence arises 
through explicit /i-dependent terms canceling against implicit //-dependence in e and 
m. In practice the scale independence only arises approximately and the goal is to 
have a procedure in which one can systematically remove scale dependence order by 
order in e. 

In the next two subsections we present simple examples to explain the standard 
renormalization group terminology from the viewpoint of light-front field theory, and 
to elucidate coupling coherence. 

3.6.1 Relevant, marginal and irrelevant terminology (on the 
light-front) illustrated through simple examples 

Because of the simple analytic form of the transformation we will use the Wegner 
transformation in this subsection. Recall Section |3.4| , especially Eq. ( 3.59|) which is 



the exact fiow equation for the effective interactions where Hg = h + Vg and 

Vsij = exp(-sAjj.) Vsij , (3.80) 

with Aij = Si — Sj. 

^^Strictly speaking, i/QED is a 'trivial' theory, the limit A — > oo should not be taken. However, 
even in the trivial scenario, as long as A <§C TOexp(|^) ^ m x 10^^°, logic prevails in all practical 
calculations. 
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From Eq. (|3.59|) we see that Vg only changes at non-hnear order and we have 

^ = Oivl) . (3.81) 

The standard renormahzation group terminology is defined with the assumption that 
the non-linear terms are small. At linear order we have 

dv<i 

— ^ = . (3.82) 

as 

Now we write a general interaction in scalar light-front field theory in 3 + 1 di- 
mensions. Since the physical idea behind renormahzation group works best for local 
interactions, and longitudinal locality is absent in a light-front approach, we only 
assume transverse locality and only keep track of transverse momenta in this ini- 
tial example (we discuss a longitudinal renormahzation group below). Given this, a 
general interaction (with no zero mode terms such as 01020304) has the form 



00 CO M—1 

Vs = 

n=0 M=2 m=l 



C„„m / levr^^^ (Pl + ^2 H Wm- Pm+l - Pm+2 Pm) 

„_n»,r_n™_i Jl,2,-,M 

xa\al ■■■al^ am+iam+2 ■ ■ ■ ^Af p1" , (3.83) 



where = Sp^= j ^16^^^+^^ ' '^^ ~ '^(Pi)? ^^e commutation relations are 



d^vj:dp+e{p+) 

[oi,o5] = 167r=^p+53(pi-p2) , (3.84) 



and the factor is a shorthand^ for a product of M transverse momenta each to 
a particular power with the sum of all the powers being equal to 2n — this is all that 
matters in order to classify these operators; note n can not be negative based on 

the assumption of transverse locality. CnmM are dimensionless couplings. We have 

^"^This is '2n' based on the kinematic rotational symmetry about the longitudinal axis with cor- 
responding generator . 
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explicitly made these couplings dimensionless by including the appropriate power 
of the cutoff through the factor s^.Q A simple exercise shows that to make CnmM 
dimensionless we must have 

M + 2n - 4 , , 

D = — —- . 3.85 

4 

Taking a derivative of this general form of the interaction givesQ 
1 dv. 



V / 16tT^S^ {pi+P2-\ hPm- Pm+1 " Pm+2 Pm) 

t t t 2n I dCnmM \ /„ 



Thus at linear order [recall Eq. (|3.82| )] the flow equation implies 



+ D c^^M = , (3.87) 

at 

where we have defined t = Ins. The range of s is to cxo which in terms of t 
corresponds to — oo to oo: running the transformation (that is integrating out high 
energy scales) corresponds to increasing t. Solving this gives 

CnmM = k exp(-D t) , (3.88) 

where k is an arbitrary constant fixed by the boundary conditions. The standard 
renormalization group terminology is 

D < — > relevant coupling , 
D = — > marginal coupling , 

D > — > irrelevant coupling . (3.89) 

'^"'^Recall the dimension of s is l/(energy)^, which as far as transverse momentum goes imphes 
s — > Sg and are the initial and final cutoffs respectively. Recall that > Sg and the "no 

cutoff limit" is s„ — > 0. 



^^Note that — 0. Recall that the free Hamiltonian h only depends on the renormalization scale 



ds 



^ and not on the running cutoff scale s: ^ = 0; and a||0) is a state of h that thus only depends on 
fi but not on the running cutoff scale s. 
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We see that at this hnear order the relevant couphngs are exponentially enhanced 
as the energy scales are integrated out whereas the marginal couplings remain fixed 
and the irrelevant couplings are exponentially suppressed — hence the terminology. Of 
course this changes when the non-linear terms are accounted for, but typically (as 
long as the couplings never get too big) this only changes the results logarithmically 
in s, that is linearly in t: the non-linear changes occur slowly. 

The results of this initial example correspond to the usual equal-time renormaliza- 
tion group terminology. Recalling Eq. (|3.85| ) and unraveling the above results we see 
that there are a finite number of structures — only taking into account the transverse 
scales recall — that are relevant and marginal whereas there are an infinite number 
of structures that are irrelevant. The only allowed relevant, marginal and irrelevant 
structures according to the transverse renormalization group (with no zero-modes) in 
scalar light- front theory in 3-|-l-dimensions are respectively 

D < — > (ctia2, a\a2a3, alala^) , 

D = — > {a[a2 p\, a\a2a^ciii, a\a2a^a4^, a\a2a\ai) , 

D > — >• (a|a2 a\a2a^ 01020304 p^, 0I02O3O4O5, . . .) . (3.90) 

Now we work the previous example, scalar light-front theory in 3 + 1-dimensions, 
but only keep track of the longitudinal scales. We will not assume longitudinal locality. 
A general interaction (with no zero mode terms such as O1O2O3O4) has the form 

00 00 M—l „ 

""s = E E E '^riniM / IGlT^S^ {pi+P2-\ ^ Pm - Pm+1 - Pm+2 Pm) 

n=-ooM=2m=l J1,2,--,M 

xa{al ■ ■ - oj^ o^+iO^+2 ■ ■ - Om (p^T > (3-91) 

using the same notation as in Eq. (|3.83|) , but with p^" — ^ (p^)" (note now —00 < 
n < 00). Note that s — > l/(energy)^ — > {p^)^, so now to make the coupling CnmM 
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dimensionless we must have 



D = ~. (3.92) 



'M' does not enter (unlike the transverse renormahzation group) because a\ are di- 
mensionless with respect to the longitudinal scales. The discussion now follows as 
below Eq. (|3.85| ) except with this new value for D. Thus we have: Linear order 



longitudinal renormahzation group implies 

CnmM = kexp{—Dt) = kexp{nt/2) , (3.93) 

where k is an arbitrary constant fixed by the boundary conditions. Recall that in- 
tegrating out energy scales corresponds to increasing t. Sticking to the standard 
terminology of Eq. (|3.89D we have 



n > — > relevant coupling , 
n = — > marginal coupling , 

n < — > irrelevant coupling . (3.94) 

We see that the sequence of nonlocal operators (^^^ <P^, 0^, • • • becomes more 

and more irrelevant in the longitudinal renormahzation group sense, whereas recall 
that in the transverse renormahzation group sense the sequence {p^^ 0^, (<9^) 0^, 
. . . becomes more and more irrelevant. The same type of simplification that occurs 
in the transverse renormahzation group when locality is assumed may occur in the 
longitudinal renormahzation group when non-locality is assumed. This may be a 
deep conclusion, with perhaps far-reaching consequences, but complications occur in 
practice when running a longitudinal renormahzation group , and it is not clear at 



present how to proceed. It is clear however, with these last two examples as witness, 
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that longitudinal and transverse scales must be treated differently in light-front field 



theory as long advertised by Wilson and collaborators |T4|. 



3.6.2 Simple example to elucidate coupling coherence 



This example is from |30 . However, we will use the G-W transformation here. 



Recall the flow equation in the G-W formalism 

v^ = v + 8v^- /^K, T^>]d\' . (3.95) 

This example is massless coupled scalar theory in 3 + 1 dimensions: 

V = J d^xHx' + + 11^0^0^) . (3.96) 

We normal order this interaction and drop zero-modes. 

To elucidate the discussion on the renormalization scale /x at the beginning of this 
section, we show how it enters the calculation through the running of the marginal 
couplings at one- loop. The counterterms for these above three marginal couplings are 
defined by 

K = - r[vx',Ty]dX' , (3.97) 

J A 

where A is the initial bare cutoff and fi is the renormalization scale. Given this, the 
above effective interaction at the final cutoff A for these three marginal couplings is 

vx = v- r[vy,Ty]dX' , (3.98) 

where v has implicit dependence on /i through the couplings gi, g2 and g^. Note that 
the bare cutoff dependence is gone. 

Now we calculate the two-particle goes to two-particle matrix elements of vx 
through one loop, setting the transverse Jacobi momenta to zero — which isolates 
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the running of the marginal couphngs — and derive how the marginal couplings gi, g2 
and gs must run so that the respective matrix elements of vx are independent of the 
renormalization scale /i. We will not go through the explicit calculation, but will just 
present the results. Through one loop we obtain 



where i = ln(//). /i has dimension mass and arose from a light-front energy renormal- 

2 

ization scale now integrating out energy scales corresponds to decreasing t. We 
see that the theory is not asymptotically free. 

Coupling coherence enters when we search for solutions to these flow equations. 
The initial theory has three independent marginal couplings. Under what conditions 
will this number be reduced? Also, the initial theory is not 0{2) symmetric with re- 
spect to 4>1 + (j)l. Is it possible to somehow restore this symmetry? Coupling coherence 
is one way. The first postulate (the second one comes below) of coupling coherence in 
this example is that g2 and g^ do not depend explicitly on the renormalization scale 
II but rather only implicitly through their dependence on gi — one can say they run 
coherently with gi. In equations this postulate is that ^ = and ^ = 0. Thus, 
through one loop we have 




(3.99) 



dgi 
dt 
dg2 
dt 



dgi ^ 
dt 

dg2 dgi 
dgi dt 




3 




(3.100) 
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where we inserted the results from Eq. ( p. 991) on the right-hand-side of these equations. 



Inserting the top equation into the other two gives 



dgi 

^^{al + gl) = 1 (9193 + 9293 + ^91) , (3.101) 

once again valid through one loop. 

Two non-trivial solutions to these two coupled equations given the boundary con- 
dition of coupling coherence {92\g-^ >o ~ ^ and g-^l^^ = 0) are (i) a decoupled 



theory: 



{92 = 91, 93 = 9i) — ' = 77 



9i \ f4>i + . (4>i- 4>2\^ 



4! [V 23 ; V 23 ) 
and (ii) an interesting 0(2) symmetric theory: 

As promised, coupling coherence is a way to reduce the number of independent cou- 
plings and restore symmetries not manifest in the flow equations. This choice of which 
marginal coupling is to be the independent coupling is where all the physics lies. In 
QED it is clear that e is the coupling of choice. 

For further simple one-loop examples which show how coupling coherence leads 
naturally to massless gauge bosons when explicit gauge invariance is broken by the 
regulator, to massless fermions when explicit chiral invarianceQ is broken, and to the 
hidden cj) — > —cf) symmetry in scalar theory when working in the broken phase (of 
particular interest for light-front calculations without zero-modes), see . 



^•^ Chiral symmetry is very interesting on the light-front. Its charge measures helicity. See 
and references within for further discussion. 

51 



CHAPTER 4 



POSITRONIUM'S GROUND STATE SPIN SPLITTING 



Now we apply the procedure outlined in Sections |3]^ and ^]3| to obtain positro- 
nium's ground state spin splitting through order a^. First, we derive Hx to second 
order in e. This includes a discussion of the effective fermion self-energy, but the effec- 
tive photon self-energy (for A < 2m) and electromagnetic coupling do not run at this 
order. Then we move on to the diagonalization of H\. This starts with a discussion 
of our zeroth order Hamiltonian, TYq, which will be treated analytically. We introduce 
a coordinate change that takes (x G [0, 1]) — > (k^ G [— cxd, cxd]), which allows easier 
identification of Tio- We solve for the spectrum of Tio exactly, which among other 
things, fixes the a-scaling of the momenta in the matrix elements in BSPT. Then we 
move on to a derivation of the perturbative effects coming from low-energy (energy 
transfer below A) photon emission, absorption and annihilation at order e^, which 
includes a discussion of the full electron and positron self-energies and a derivation 
to order of the complete exchange and annihilation interactions. Given this, we 
determine the range of A that allows the effects of low-energy (energy transfer below 
A) photon emission and absorption to be transferred to the effective interactions in 
the |ee) sector alone, and at the same time, does not cut into the nonperturbative 
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features of the solutions of T-Cq. Finally, we proceed with BSPT in V noting that all 
shifts appear in the few-body sector, |ee), alone. 

We present the well-known result from equal-time calculations to compare with 
our results later. Through order a'^ an energy level in positronium with quantum 
numbers (n, /, S, J)Q according to QED is given by 



EnlSJ - + 



where Ryd = ™^ ~ 13.6 eV 



11 / 1 \ 1 



32n4 V 21 + lJ V? 



o? Ryd 



_ 7 1 - (5z,o / 3/ + 4 1 , 3[-l \ 

ei^j - -di,o + ^^^^ l,^^-'+^(/ + l)(2/ + 3) " ^''^WTT) ~ ^'''-'J{2r^) ) 

and Eioj = 0. This is a well established result first derived in 1947 0. For the ground 
state spin splitting this gives 

E,^s -E,^s = I Ryd (4.1) 



through order in the standard spectroscopic notation, n^'^^H. In passing it is 
interesting to note that in positronium (since fine and hyperfine structure are at 
the same order) no degeneracy with respect to J remains through order — unlike 
hydrogen where through order there is the famous 25"! and 2Pi degeneracy [and 
it is of course this splitting at order a^ln(l/a) that is the famous Lamb shift, the 
dominant part of which is calculated in Chapter Q . 



^^n is the principal and I is the orbital angular momentum quantum number; S is the total electron 
plus positron spin quantum number and J — 1+S with J the corresponding quantum number. 
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4.1 Derivation of Hx through second order 



From Eqs. ( p.l5| )-( p.48| ), the final renormalized Hamiltonian to second order is 
given by 



(i\Hx\j) = fxij I hij + Vij + VikVkj 



9ikj ^ 

^ik 
'jki 



+ ^ 1 + Sv^ + 0{e') \ . (4.2) 



gl^^^ is written in Eq. ( p.47| ) and v is given by Eq. ( p.51| ) with all the proton fields 
set to zero. 

4.1.1 RenormaUzation issues in positronium 

The form of Sv^ follows from the constraint that Hx satisfies coupling coherence. 
To order the fermion and photon masses run, but the coupling does not. First, we 
discuss the result for the electron self-energy coming from the second-order effective 
interactions in Hx- Consult Appendix |B| for a listing of the free matrix elements of 
the canonical interaction v used in this dissertation. Specifically, we calculate a free 
matrix element of Hx given in Eq. ( [4.2|) in the electron self-energy channel. 

At second-order, an electron self-energy effective interaction arises because photon 
emission in Hx is restricted to be below A. That is, the energy scales A down to A 
from photon emission have been 'integrated out' and placed in effective interactions. 
At second-order the explicit form of this electron self-energy effective interaction is 
(note Ei = Ej so fxij = 1) 

e-^(2) , ., (2), 



electron self— energy 



Pi (Aj) 

E / ^(Pt - (^'P^)0{p^ - eV^){p3S3\Vee-y\PeSeP'rSy){peSeP^S^\Vee^\piSl 
\o\i) Jpe.P^ 
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X , K 



P^s^ 



1 -X , - k: 



Figure 4.1: Momenta and spin label conventions for the second-order effective electron 
self-energy interaction. 



X 



167T-^6\pe +P^- Ps) IQTT^S^Pe + - Pi) {f Mk " f\jk) I ^jk + 



r (2) 

+ 1 

Pi 



where Vee-y = / (fx^dx TCee-y is the canonical emission and absorption interaction, 
and the initial, intermediate and final free states are labeled \j) = |1,2), \k) = 



\PeSe,P'yS-y,P2S2) and \i) = |3,4) respectively. See Figure ^TT] for the momenta and 

(2) 

spin label conventions. 6vj^ are the aforementioned second-order counterterms to be 
determined through coupling coherence below, 'e' is the infrared regulator discussed 
in the paragraph containing Eq. ( p.67|) that we are forced to introduce. We define 
our Jacobi variables by 



V = P1+P2 , 
Pi = {xV^ , K + xV 

Pe — {xePi 1 l^e + X^Pi 
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1 - Xe)pt , - + (1 - Xe)pi 



Note that 



Pe +P~y 



Pi' + + li) 



Pi 



(4.3) 



The above becomes (see Appendix ^ for the matrbc elements of Vee^ in the free basis) 

2x. 



(2) 



rl — e/x 



e/x 167r'^a;e(l — Xe) 

{f\jk - fxjk) «e + ^ + 2;e(l-Xe) 



2 A I m2(l-xe)2 



1—Xe 



(4.4) 



where 



fx'jk = d [ A'^x + - 



•^e(l •^e) ^ 



with 



GA'[Xe] = Xe(l — Xe) ( A'^X + 



Including these constraints we havef^ 



5S 



(2) 



^2 ^a;^ rf(Kj 



2+A2a; 

2 / 2 



1-e/x 



+ ^ + 



X e ( 1 3^' e ) 



m2 -f A2a: 

I m^(l-a:e)^ 



(i(K 



K-l+m^ _|_ Kj _ ^2 



+ 5v 



(2) 



Performing this integral gives 



(4.6) 



2x, 



IQ'K^Xei'i — Xe) 

(4.7) 



^xv(^) I X 



(4i 



^^We use e/x <^ 2^\,'i ■ We see that a non-zero renormalized electron mass m "regulates" the 
infrared electron momentum [pf — > 0] divergence. 
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where we have defined 



.^(2) _ a 
Ztt 



m + A X 



X"x 



(4.9) 



For respective ST^y terms of Eq. (|4.8| ), A' = A and A. Note that the energy dependence 
on the electron's relative transverse Jacobi momentum, k, does not change. Also, note 
that the result is infrared singular. 

Now we determine 5vf^ through coupling coherence given the above results. Con- 
straining the electron mass to run coherently with the cutoff according to Eq. ( p.76| ) 
amounts to the requirement 



- 5Sf ) + 5^'' + O (e^)] = [m^ + 5v'"^ + O (e^) 



This fixes the countertermQ 



5v 



(2) 



(4.10) 



(4.11) 



and to second order the fermion mass renormalization is complete. 

In summary, through second-order the coherent electron mass-squared is 



ml = 



m — (oL^ — ol^x ) + "^"^A 



(2) 

A 5 



(4.12) 



where SH^^ is written in Eq. ( [4.9|) . In H\ there is of course still the photon emis- 
sion interaction below A that must be considered. The form of this interaction is 

fx J (fx-^dx'Hee-y This is considered below in Subsection 4.2.2, and the resulting 

36 Any finite ©(e^) running cutoff independent term can be added to the counterterm, and 
Eq. (4.1C) would still be satisfied. This term can only depend on the renormalization scale /i. 



Setting this term to zero is our choice of renormalization prescription. 
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combination of m\ and these effects from "low-energy emission" add to the physi- 
cal electron mass-squared fn^i^y^^ which as will be shown is infrared finite and with 
our choice of renormalization prescription (see footnote through second order, is 
given by m^, the renormalized electron mass-squared in the free Hamiltonian h.^'^ 

For arbitrary A, the photon mass also runs at order e^. The discussion follows that 
of the electron mass except for the fact that the running photon mass is infrared finite. 
For < (2m) ^, the resulting coherent photon mass vanishes because pair production 
is no longer possible. In this chapter we choose A^ ^ m^, thus the photon mass 
is zero to all orders in perturbation theory, as required by gauge invariance. There 
are additional difficulties with dependent marginal couplings that are encountered at 
0{e^), but this is beyond the focus of this dissertation. 

4.1.2 Hx through order e^: exchange and annihilation chan- 
nels 

To complete the derivation of H\ through second-order, we need to write the 
coherent interactions for the exchange and annihilation channels in the |ee) sector. 
At second order, these come from tree level diagrams, with no divergences or running 
couplings. Thus, the coherent results follow from 

6vT = - r[vy,T^>]d\' . (4.13) 



(^) _________ 

To show that 5vp^ produces a coherent interaction recall Eq. ( ^.341 ). We have 



= V- r[v';\Tl?]d\' + 0{e') (4.14) 

J A 

^^So this is a "physical" renormahzation prescription; not that it matters, because physical results 
are of course independent of the choice of renormalization prescription. 
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Figure 4.2: This illustrates the spin and momenta label conventions used for this 
positronium calculation. 



- r[v';:,Ti']dX' + Oie')\ , (4.15) 

which satisfies the coupling coherence constraint, Eq. (|3.76| ). At second order this 
seems trivial, but at higher orders the constraint that only e and m run independently 
with the cutoff places severe constraints on the Hamiltonian. 

Given this second-order interaction, the free matrix elements of H\, shown in 
Eq. ( [4.2| ), in the exchange and annihilation channels are 
Exchange Channel: 

(e(3)e(4)|i/,|e(l)e(2))|_,_ 



X, exchange 



1Qtt'^5'^{Pi + P2 — P3 — P4)yxx'{l — x)(l — x') 



(4.16) 



where 



Ml 



X 



'^^(|Ai|-|A2|)^^(|Ai|-A2) , 0(|A2|-|Ai|)0(|A2|-A2) 



DENi 



+ 



DEN2 



(4.17) 
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V2 = -e'e{>^- 



Ml-M',' 



. L/.tJi fin 



The variables are defined as follows (see Figure [4.2| also): 

(xV^, K + xV^) , P2 = ((1 - -K + (1 - x)V^) 

[x'V+, k' + x'V^) , P4 = ((1 - x')V^, -k' + (1 - x')V^ 

(x — x'Y 



(4.18) 



Pi 

P3 



Si 



xx'(\ — x)(l — x') 

S2 



A- A 



X — X' 



1 — x' 1 — X 



5 -^2 



X — X' 



X 



DENi 



K-' + i s ejk k!" (s = ± 1 and s = -s) ; eu = -t2i = 1 , en = 622 = 



A2 



XX' 



XX' 



KyY — 

x(l -x) ' ° x'(l -x') 



Annihilation Channel: 



(e(3)e(4)|/f,|e(l)e(2))| 



annihilation 



X,annihil 



167r'^(5^(pi + P2 — P3 — Pi)\lxx'{l — x)(l — x') 



(4.19) 



where 



e (mi - M'o^) e{Ml-x 



Ml 

M''-Mt]9(M'r:'-X 



f2 



Ml - M',' 



) ^SrS2^. 



1 



(4.20) 
(4.21) 



and 

N2 



'^SlS2 '^5453-^3 ' -^4 ~l~ ^SlS2'^SlS4'^S3S4 



2m2 



xx'(l — x)(l — x') 
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^ 1 — x' x''^ 1 — X X 

V2 and V4 are canonical instantaneous exchange and annihilation interactions, respec- 
tively, with widths restricted by the regulating function, fx. Vi and V3 are effective 
interactions that arise because photon emission and annihilation have vertices with 
widths restricted by the regulating function, fx- 

4.2 Diagonalization of Hx 

First we discuss the lowest order spectrum of Hx, after which we discuss BSPT, 
renormalization and a limiting procedure which allows the effects of low-energy (en- 
ergy transfer below A) emission to be transferred to the |ee) sector alone. 

4.2.1 7^0, a coordinate change and its exact spectrum 

Tio, as motivated from the form of our second-order effective Hamiltonian, in the 
|ee) sector is 

Ho = h + V^ , (4.22) 



where h is the free Hamiltonian given in Eq. ( |3.11|) , and V^^ is given by [using the 



same variables defined below Eq. (|4.18|) ; note, Kz is defined below by Eq. ( [4. 241) ] 

(e(3)e(4)|V,|e(l)e(2)) 



1Qtt^5^{Pi + P2 — P3 — Pi)yxx'{l — a;)(l — x' 



where V, ^ - . (4.23) 

In all other sectors, we choose T-Cq = h. This is convenient for the leading order 
spin-splitting calculation (order a'^) in positronium, but for example in the Lamb 
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shift calculation it is more convenient to keep the Coulomb interaction between the 
oppositely charged valence particles to all orders in all sectors. 

As already mentioned, Tio in the |ee) sector is motivated from the form of our 
second order effective Hamiltonian Hx: it arises as the leading order term in a non- 
relativistic limit of the instantaneous photon exchange interaction combined with the 
two time-orderings of the dynamical photon exchange interaction. We choose it to 
simplify positronium bound-state calculations. Other choices are possible, and must 
be used to study problems such as photon emission. Later, in BSPT this choice 
is shown to produce the leading order contribution to positronium's binding-energy 
(order a^) as long as ma^ ^ A ^ ma. 

The coordinates Hz and in Eq. ( [4.23|) follow from a standard coordinate trans- 



formation that takes the range of longitudinal momentum fraction, x G [0, 1] to 
Kz G [—00,00]. This coordinate change is 

^ = 1+ (4.24) 

2 2^/ k"^ + k1 + m? 



We introduce a new three- vector 

P = (/«,/«.)• (4.25) 

Note that 

Ml ^ 4^ = 4K + P') (4.26) 
x[\ — X) 

is invariant with respect to rotations in the space of vectors p. The nonrelativistic 
kinematics of Eqs. (|3.2|) and ( p.3|) in terms of this three- vector become 

^ = 0(a) . (4.27) 
m 
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Note the simple forms that our "exchange channel denominators" take in the 
nonrelativistic limit 



DEN, = (p-p-r-("--"-^(P'-P^'^+0 



DEN2 = (p-pO + 



m 



m 



+ 



P \ 2 
— I m 

m 

P \ 2 

m 

m 



(4.28) 
(4.29) 



Also note the form that the longitudinal momentum fraction transferred between the 
electron and positron takes 



, K, - 4 (p^ 4 - p^/t,) 

X — X = — z 1 : — 5 h O 



2m 



(4.30) 



These formulas are consistently used throughout this chapter. 

Now we describe the leading Schrodinger equation. We seek solutions of the 
following eigenvalue equation 



(4.31) 



where S 



N 



V+ 



Hq is diagonal with respect to the different particle sectors, 



thus we can solve Eq. ( U-31\) sector by sector. In all sectors other than |ee), Tio = h, 
and the solution is trivial. For the |ee) sector, a general I'lpNi'P)) is 



X (t)N{xKSiS2) bl^ipi) dl^{p2) |0) 



(4.32) 



with norm 



{^NiV)\ijN'iV')) = 5nn'IQ7t'V^S' (V - V) 
Jd'^nJ^dx- 

=^ E 77"^ (f)*]sf{xKSiS2)(j)N'{xKSiS2) 



S1S2 



6nn' 
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The tilde on (p^ will be notationally convenient below. In the |ee) sector, Eq. ( 4.31 ) 
becomes 

[mi - 1^^) Mx'n's.s,) = E I^^^ V, Mx.s,S2) . (4.33) 

After the above coordinate change, this becomes [J(p) is the Jacobian of the coordi- 
nate change written below] 



SlS2 



d^pp{p)J{p') 

167r3 



0Ar(pSiS2) 



(4.34) 



where the tilde on the wavefunction has been removed by redefining the norm in a 
convenient fashion 



5mn' 



SlS2 



SlS2 



= J2 d^P 4'*NiPSlS2)(j)N'{pSlS2) , 

SlS2 

and the Jacobian of the coordinate change is 



(4.35) 



dn, 2(p2 + m2) 
Note that the Jacobian factor in Eq. ( [4.34|) satisfies 



(4.36) 



J{p)J{p') 



1 
2m 



P' 



2nl + p'^ + 2«:f 



Am? 



+ 



' 4 /4 



(4.37) 



Before defining TYq in the |ee) sector we mention a subtle but important point in 
the definition of TYq- "^o ^he |ee) sector will not be defined by Eq. ( U-34 ). Rather, it 
will be defined by taking the leading order nonrelativistic expansion of the Jacobian 
factor in Eq. (|4.34|) . This gives 



M'^ - 4(m^ + p 



/2^ 



^n{p's3S,) = E / ,^'3"^^ V; MPS1S2) , (4.38) 
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where is defined in Eq. ( |4.23| ). This Ho will be diagonalized exactly, and the 



subsequent BSPT will be set up as an expansion in V = H\ — TYq, which will then be 
regrouped in terms of a consistent expansion in a and Q;ln(l/a) to some prescribed 
order. First, we discuss the exact diagonalization of Hq. 

Putting the expression for into Eq. ( ^.381 ) results in the following equation 

(~B^ + ^] Mp'ssS,) = ^ I ^^-^Mps^s,) . (4.39) 



m J l-K^ J (p — p'j^ 

This is recognized as the familiar nonrelativistic Schrodinger equation for positronium. 
Note that we have defined a leading order binding-energy, —B^^ as 

M\ = 4m2 + 4mS^ , (4.40) 

where is the leading order mass-squared. Note the difference in the definition 
of our leading order binding-energy and our exact binding-energy as given by = 
(2m + -Bjv)^ [see Appendix ^ for further discussion on this difference], where is 
the exact mass-squared. 

To proceed with the solution^ of Eq. ( 4.3S ) note that there is no spin dependence 



in the operator so the spin part just factors out 

0i/,^e,s_(p's3S4) = 0z.(p')^SeS3^s_^4 • (4-41) 

We rewrote as [y^ Se, s_), where (se, s_) label the spin quantum numbers and u labels 
all other quantum numbers, which are discrete for the bound states and continuous 
for the scattering states. 
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The solution of Eq. ( |4.39| ) is of course well-known, but introducing hypersphcrical harmonics 



(which may not be so well-known) is essential for the later analytic calculation in second order 
BSPT, so we go through some detail here. 
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The solutions to Eq. ( 4.39 ) are well known. For Bj^ < 0, following Fock |39|, we 



change coordinates according to 



mB^ ^ -el, (4.42) 
u = (mo, u) , (4.43) 

uo ^ cos(^) = 4^' (4-44) 

u = — sin(to') = sin(u;) [sin(6') cos(0), sin{6) sin(0), cos 
P 



2e„p 
e2 +p2 



(4.45) 



Useful relations implied by this coordinate change are in Appendix 0. Note that in 
our notation we anticipate that p will be given by {n,l,m^), the usual principal and 
angular momentum quantum numbers, and that the leading order binding-energy will 
depend only on the principal quantum number, n. Given this, Eq. ( [4.39] ) becomes 

M^P') = / j:7^M^p) , (4.46) 



n 



\u — w 



where 



/ 2 2 \ 2 

4(enJ 



2 



Using Eq. ( p.20| ) of Appendix 0, this is seen to have the following solution: 

M^,)=Y,{np) with — = 1, (4.48) 

where Yj^(flp) is a hyperspherical harmonic. Thus, 

ma ^ ma^ 



" 2n " 4n2 

and 



= > (4-49) 



4 (en 



(P) = ..2 ':2^2 ^^(^p) • (4-50) 
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This is the standard nonrelativistic solution for the bound states of positronium. 
This completes the solution of TIq for the bound-states. The scattering states are also 
needed in our second-order BSPT calculation. We use Green's function techniques to 
include the scattering states where required. 

4.2.2 BSPT, renormalization and a limit 

Here we use the BSPT formulas (appropriately generalized to the degenerate case) 



of Section to analyze positronium's leading ground state spin splitting. The 



potential to be used in BSPT is 

V = H^-H,, (4.51) 
where the eigenvalue equation for TYq is given by Eq. (f4.38|) , and H\ to second order 



is given in Section [4.1| . We will be perturbing about the nonperturbative eigenstates 
of 7^o• 

First, we discuss electron mass renormalization. In second-order BSPT there is an 



electron mass shift coming from the f\v part of Hx, with v given by / (Px-^dx Ti, 
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[see Eq. ( |2.53| )]. This is photon emission and absorption restricted by the regulating 
function, fx- The calculation is similar to that of Subsection ^4.1.1| . Assuming {M.\ — 
M-D = O(e^), this electron mass-squared shift is 

6m^ = -(5Sf + O(e^) . (4.52) 

6T,x is the same function that was defined in Eq. (WM). Checking the consistency: 



using this result [Eq. ([4.52|) ] one obtains (A^^— A^q) = (4m^+4mi3^ — 4(m^ + p^)) = 
(9(e^), and our initial assumption is satisfied. Combining this with the coherent 
electron mass-squared m\ of Eq. ( |4.12|) , the physical electron mass-squared m^^^^ 
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through second-order is 

mlhys = ml + Sm'^ 

= - {6J:'1' - 5Sf ) + 6vT] + [-^Sf] + 0{e^) 

= + O(e^) . (4.53) 

In this last step we recalled the result from coupling coherence of Eq. ( [4.11j ). We 
see that through second order, the physical electron mass is equivalent to the renor- 
malized electron mass in the free Hamiltonian h. As already mentioned, we see that 
our choice of renormalization prescription mentioned in footnote ^ corresponds to a 
physical prescription as the physical electron mass is equivalent to the renormalized 
electron mass through second order. Of course we do not have to use this physical 
prescription, but it is our choice here, and must be maintained consistently in higher 
order calculations. As promised below Eq. (|4.12| ), we see that treating photons per- 
turbatively has led to an exact cancellation of the infrared divergence in the coherent 
electron mass-squared m\, and the physical electron mass through second-order is 
infrared finite. 

Now we move on to the discussion of BSPT. The only channels to order are 
exchange and annihilation. Parts of these effective interactions are given in Subsec- 



tion [4.1.2| . We also need to include the perturbative mixing of the \ee'y) and I7) 
sectors with the |ee) sector arising from fxv, with v = J (Px-^dx'Tiee-y- In second- 
order BSPT this gives rise to the following effective interactions that must be added 
to K , and K ^ , of Eqs. (14.161) and (|4.19|) respectively. 

A, exchange \,anmniL ^ \| v ^ 



Exchange Channel: 
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with 



DEN^ 



[k - K'f + l{x- x')A +\x- x'\ (1 [Ml + M',') - Ml) 



and 



A 



1 — X 1 — x' 



X' 



Annihilation Channel: 



M^. 



(4.55) 



Note that in a nonrelativistic expansion [after the coordinate change of Eq. (|4.24| )j 
the above "exchange channel denominator" becomes 



DEN, = (p - pO' + \x- x'\ {Ml + M',') -Ml)+0 



2 

— m 
mj 



The full exchange and annihilation channel interactions to order are 



(4.56) 



V 



exchange 



X, exchange ^''^ ' 



'^annihil '^X.annihil ^''S 



(4.57) 
(4.58) 



where Eqs. (14. 161) and ([4.191) give V, . and V, . , respectively. 

1 \| 1/ \| V \, exchange X,anmhii ^ ^ 

One way to summarize the results, recalling the form of Eq. (|4.34| ) and the norm in 
Eq. (|4.35|) , is to state: The full order effective interactions give rise to the following 
first order BSPT shift of the bound-state mass-squared spectrum of I-Lq 



(5*''M^(s3,S4;Si,S2) 



\y\(^n,l,mi,sx,s^ 



J £p C/V0n,i,mi(p')^(p', P, Si, S2)0n,i,mi (P) ) (4-59) 
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J Mm (v +v ) 

V IOTT"^ IGTT"^ ^ exchange anmhil J 



where 

l^(p', S3,S4;p, Si,S2) 

V lU/l " lU/l " ^ 

The Dirac notation in Eq. ( |4.59| ) will be used in the remainder of this chapter. See 
Eqs. ( [4. 231 ), ( [4.57D and ( [4. 58] ) for V^, V.^^^^^^^ and V^„„i^i, respectively. The interaction 



V must be diagonalized in the degenerate spin space following the standard rules of 
degenerate BSPT. Note that in order to get all the shifts through a consistent order 
in a (in this case through order a^), V needs to be considered in second order BSPT 
also. 

The diagonalization of V in the degenerate spin space follows shortly, but first we 
determine the range of A that allows the effects of low-energy (energy transfer below 
A) photon emission and absorption to be transferred to the effective interactions in 
the |ee) sector alone, and at the same time does not remove the nonperturbative 
bound-state physics of interest. This range is 

\Ml-i2mf\ A2 

, 4.61 



'P^ ^photon 

where is the bound-state mass-squared and is the dominant energy of an 

emitted photon. After the solutions of Tio are known the a-scaling in all BSPT matrix 



elements is known and the bounds in Eq. ( |4.61|) become 



m^a^ < A^ < m^a . (4.62) 

This is satisfied under the following limit 

A^ — > a fixed number , (4.63) 
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rr?a 



A2 



oo 



(4.64) 
(4.65) 



Given the nonrelativistic limit 



m 
A2 



oo 



(4.66) 
(4.67) 



this imphes 



(4.68) 



with 



2 < < 4 



(4.69) 



Note that this "window of opportunity" is available to us because, (i) we have an 
adjustable effective cutoff ^ in the theory, and (ii) QED is a theory with two dy- 

2 2 2 

namical energy scales, and ^^^p, a fact known for a long time, and the reason 
that QED calculations have been so successful over the years. 



Given the above limit [Eqs. ( [4.63|) -( ^.65|) ] 



Q (}? - 4|p2 - p'^l) , Q (I All - A^) , (I A^l - A^) ^ 1 
e (4(p2 + m") -X^),e (4(p'' + m^) - A^ 
d[x^-\Ai\),0[x^-\A2\) —^0 
9 (a' - 4(p2 + m^)) , (A^ - 4(p'' + m^) 







(4.70) 
(4.71) 
(4.72) 
(4.73) 



Now we proceed with the diagonalization of V in the degenerate spin space [see 
Eqs. ( |4.59| ) and ( [4.60| )]. In BSPT with V we will calculate all corrections through 
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order that give rise to a spin splitting structure in the ground state of T^o- First, 
we write the general V more explicitly given the above limits in Eqs. ( [4.70| )-( |4.73| ): 

^(P', -SS, S4; p, Si,S2) 



IGtt^ 2m 



+ 2kI + p'^ + 2Kf ^ ^ / p^^ 



X 



^ ^ -^c , (4.74) 



IGtt^ 2m 



where 



^ - ^12 ^ ^21 ^ 0^^ = [DENi - DEN2) (4.75) 



DENi DENi DEN2 

(4.76) 



e [Ml - M'o^) e [M^^ - Mo'' 



DEN^ Ml M'o' 

Note that we have expanded out the Jacobian factors according to Eq. ( ^.37| ). Also, 
DENi and DEN2 are defined below Eq. ( |4.18| ) and written in their expanded version 
in Eqs. ( [4. 281 ) and ( [4.29] ) respectively. Finally, A^^i and N2 are written below Eqs. ( [4. 181 ) 
and (|4.21| ) respectively. 

Since the eigenstate wavefunctions of Tio force p to scale as p ~ ma, it is useful to 
note the a-scaling of the matrix elements of V in momentum space in a nonrelativistic 
expansion. Recalling that we are always assuming a — ^ (without which our matrix 
elements would not have a well-defined a scaling), we see the following structure 
arising 

V = V'"'' + V^"^ + V'"' + ■■■ , (4.77) 

where a momentum space matrix element of V^*^' scales as V^^^^ . Thus in first- 
order BSPT these respective terms scale as (this follows noting / d^p\(pN{p)f = 1 
and V is not diagonal in momentum space) 

S'^'mU' = (07v|v"ViV') ~ a^-"' . (4.78) 
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To be consistent through order in first-order BSPT we need to look at all the matrix 



elements of V^'^' with S < l;0and in second-order BSPT, since lAif — M.n\ 



2 2 

~ ma, 



we need to consider V"''^' in second-order BSPT (since a^a^/a 



3^,3 /^,2 



a 



A final discussion that we must have, before we write out these expressions for 



V , is how we are going to deal with DEN4 and DEN^ defined above.Qj These 
denominators are dealt with by noting the following formulas 



e{a - h) ^ e(h-a) 



+ 



i e{a-h) + e{h 

2 a 

ie{a-h)-e{h-a) I e{a - h) - e{h ~ a) 
2 a 2 h 



_a) i e{a-h) + e{h-a) 
h 



1/1 1\ 1 \a-h\ n 1 
2 Va ^ 6/ ^ 2 a-h U ~ & 

-(- - 



1 |a-6| 



2 ah 



(4.79) 



To proceed it is useful to note 
DEN^ = (p - - 



(/s:^-0(p2-p 



/2\ 



m 



DENi = DEN2 



1 1 

+ 



2(/s:,-<)(p2-p 



+ 



o 



m 



m 



2 \DENi DEN2 



m 
1 



m 



m 



(4.80) 
(4.81) 



+ 



P-P 

^(«:.-<)2(p2-p'^)2 



+ 



P 



.(4.82) 



m^(p — p')^ 

Especially note that this last equation scales as: + 1 + a"^ + ■ ■ ■, i.e. the corrections 
start at order 1 (not order ^); this implies that only the ^p_p,-^2 term of Eq. ( [4.82|) 



39 



e P 



For example, i^p^p,yi 



5 = 



Actually the DEN5 term is handled with analogous techniques as the DEN4 term, and has 
even smaller corrections than those of DEN4. Thus, we will just discuss the DEN4 term in what 
follows and here state the result for the DEN5 term: Take DEN5 — > 4m^; the corrections to this 
start shifting the bound state mass at order a^. 
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contributes to the spin splittings to order a^. But we still have to discuss the second 
term that arises in Eq. ( [4.79| ). This term is given by 



DEN^ 



llDENi- DEN2\ 1 



2(k,-4)(p^-p'^) 

m 



2 DEN1DEN2 2 fp-p')^ 

second term ^ ^ \tr tr / 

Including A^^i, this starts out as an 0{q^) spin conserving contribution, which does not 
contribute to the splitting. The next order contribution is 0{a^) with spin structure, 
but is odd under p ^ p', and thus integrates to zero in first-order BSPT (in second- 
order BSPT it contributes a^a^/a^ ~ a^). However, the 0{a^) spin conserving term 
appears to lead to an order shift to the spin splittings in second-order BSPT when 
the cross terms with V^^°' of Eq. (f4.85|) are considered;^ however, these cross term 



+0(oP) . (4.83) 



contributions add to zero due to the facts that the 0{a^) term including Eq. ( [4.83|) 



conserves spin and the 0{a^) term including Eq. ( |4.83|) is even while the term from 
Eq. ( [4.85|) is odd under p ^ p'. 



To summarize the preceding discussion of DEN^ and DEN^, we can say that 
through order ct^, for the spin splittings of positronium, there are no relativistic 
corrections to the following replacements: 

DEN^ — >{p- p'f and DEN^ — > Am^ . (4.84) 

This is valid for the ground and excited states, but in what follows we specialize to 
the ground state for simplicity. 

Given this general conclusion about DEN4 and DEN^, we list the pieces of V that 
contribute to positronium's ground state spin splittings through order a^. Explicitly, 



'^^V of Eq. ( 4.85 ) comes from the first term on the right-hand-side of Eq. ( 4.82 ) combined with 
the complete next to leading order term of A^i . 
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as far as the a-scaling goes, we need to consider 



.(0)/ ./ N -Ce.e 



2 



r '(PV4;PS1S2) = , 3, ,.^ V^"KP'S3S^;PSIS2) , (4.85) 

47r'^(p — p'j^ 

where 



V 



/l(p's3S4; PS1S2) = Si{Ky - Ky) - i{K^ - k'J , (4.87) 

/2(P'S3S4; PS1S2) = S4{Ky - Ky) + i{Kx - k'^) . (4.88) 

Recall that Si = ±1 (i = 1,2,3,4) only. The only other interaction we need to 
consider is 



V (PS3S4;PS1S2) = 



o2 



4m7r^ 



^an^SlS2^SlS4^S3S4 ~l~ C-^x^S2S4,^S2Sl^S3Sl ~l~ 

(4.89) 



+ I '^on2 '^"^ (p — p')2 ] ^SlS20s3S4 



The constants c^^ and were introduced only to distinguish the terms that arise 
from the 'exchange' and 'annihilation' channels respectively, and c^^ = c„„ = 1 (as 
given by the theory) will be used in the remainder of this calculation. 

Two simplifications were made in deriving v''^\ First, we did not include terms 
that are a constant along the diagonal in spin space, because these do not contribute 
to the spin splittings. Second, we noted that terms of the following type integrate to 



zero 



(01,0,0.3.4 1 (p_p.)2 ^101,0,0.1.,) = (0,0,0) , (4.90) 

and thus were not included in the definition of V^^\ 

The ground state spin splitting through order contains contributions from 
V''^^ in first-order BSPT (v''°^ vanishes in first-order BSPT) and V^°^ in second-order 
BSPT. We begin with the first-order BSPT calculation. 
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First-Order BSPT: 



The lowest order wavefunctions are given near the end of Subsection [4.2.1| (see 
Appendix ^ for the hyperspherical harmonics). V^'^' in first-order BSPT contributes 
the following to positronium's ground state mass-squared 



5M^ 



5*"M^(s3,S4;Si,S2) 



^ / d'vd'v'j^^-^,j^^^v''\v\s,,s,-v.s,,s,) , (4.91) 



where 



N 



8ef 



and ei 



ma 



2n 



n=l 



(4.92) 



Using the rotational symmetry of the integrand we can make the substitution 



(P - P') 



'^2 



IP - P 



2 
3 



(4.93) 



After this, the remaining integrals are trivial and the splittings that arise from 
diagonalization of the SM"^ matrix in spin space are 



5Ml 



6M' 



5M: 



5M: 



2 4 

-m a 



2 4 

-m a 



2 4 

m a 



2 4 

m a 



(4.94) 
(4.95) 
(4.96) 
(4.97) 



where 



|1) 



+ -)-!-+) 



|2) 



V2 

I + -) + !- +) 



|3) 



|4) = | + +) 



(4.98) 
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|3> |4> 



1 

3 



2> 



5 
2 m^a^ 



J_ 
2 



Figure 4.3: dM"^ is the part of the ground state spin sphttings from first-order BSPT 
given by Eq. ( [4.91D . m is the electron mass and a is the fine structure constant. The 
state labels 1, 2, 3 and 4 are explained in Eq. (|4.98| ). The two upper most levels 
should coincide in a rotationally invariant theory; and after including second-order 
BSPT, they do. 



Figure [4.3| shows these results, which taken alone do not produce the degeneracies 
required by rotational invariance. 

Now we perform the second-order BSPT calculation. 
Second-Order BSPT : 

V^°^ gives rise to the following contribution to positronium's ground state mass 
squared in second-order BSPT 



5Ml = 6'"'m\ss,s,;s,,S2) 
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E E 



C \(Plfl,0,S3,S4, 



w 



(0) I 



,0,0,81,52) 



-1/^(1,0,0) 



Ml - 



(4.99) 



Recall that z/ = (n, I, mi), the usual principal and angular momentum quantum num- 
bers of nonrelativistic positronium (the "c" on the sum emphasizes the fact that the 
continuum states must be included also). The calculation of is tedious, but can 
be done analytically. This calculation is performed in Subsection 4.2.tj . The result is 
[see Eq. { ^JT^ ] 

-(3(71 + (72), (4.100) 



but 



24 



where g\ and g2 are given in Eqs. (|4.134|) and ( [4.135| ) respectively. 

Now we combine the bM"^ and ^M'^ matrices and diagonalize the result. The 
combined matrix is given by 



+ bMl 



The eigenvalues are 



-^(3^71+^72). 



(4.101) 



ml + ml 

ml + ml 

m'^ + m'^ 
1 2 

ml + ml 



' 2 4 

-m a 



2 

^24 

-m a 
3 

2 

^24 

-m a 



2 4 



(4.102) 
(4.103) 
(4.104) 
(4.105) 



and the corresponding eigenvectors are the same as in Eq. ( [4. 981 ). Figure |4.4| displays 
these results. 

These results translate to the well known result as detailed in Appendix |^: 

Btripiet- B singlet = -a^Ryd + 0(mct') . (4.106) 
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2 



1 

3 



|3>, |4> 



2> 



2>, |3> |4> 



1 

2 



5 
6 



1 > 



5 M? 



S M f + S M I 
2 m^cx^ 



Figure 4.4: The combined ground state spin splitting from first- and second-order 
BSPT in positronium through order a'^ is illustrated using the same notation as in 
Figure [4 .31 . SM"^ is given by Eq. ( [4.99|) and is calculated in Subsection [4.2.3| . The 
final combined result (on the right) corresponds to a rotationally invariant theory. 



We see rotational invariance in the degeneracies of the ground state n = 1 triplet 
levels exactly maintained through order a^. 

4.2.3 Calculation of 



In this subsection we perform the following sum analytically 



E E 

Se,Se 1/^(1,0,0) 



\V 



(0) I 



,0,0,si,S2) 



Ml -Ml 



. (4.107) 
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Recall that for the bound states, v = {n,l,mi) the usual principal and angular mo- 
mentum quantum numbers of nonrelativistic positronium. We must also include the 
scattering states of course. We do this with Green's function techniques as explained 
below. Recall that the spin factored completely out of our lowest-order Schrodinger 
equation, so the following notation is useful 

\<Pu,s.,s-) = 10.) ® \SeSe) , (4.108) 

Se,S-,U 

= E /^^P IP)(Pl ® \SeSe){SeSe\ ■ (4.109) 

To proceed, define the following Green's function for arbitrary E 

Ge ^ ^c\^MA. a 110) 

The factor will turn out to be useful. This Green's function satisfies the familiar 
Coulomb Green's function equation 

i'(p-p') = (£-5^)G,(p',p) + ^/dV^^P^, (4.111) 

where 



{p'\GE\p) = GEip',p) , (4.112) 
E - Am? 

and E = . 4.113 

4m 

Hostler and Schwinger independently obtained the solution for this Coulomb Green's 



function in 1964 We find Schwinger's form useful; the equation he solves is 
exactly the above equation with the following shifts in notation 

(Ze'^) — ^ a , (4.114) 

V / Schwinger 

m 

fnschwinger ^ — , (4.115) 

E Schwinger ^ E . (4.116) 
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His result is amended because the sum we need has E = M.\ = Arn^ —Aef (ei = ma/ 2 
recall) and does not include z/ = (1,0,0), as dictated by the usual rules of second- 
order BSPT. This subtraction of the z/ = (1, 0, 0) term amounts to the term "— ^" in 
Gin below. The details of how this arises can be seen in Eqs. ( [4.157| )-( [4.159| ) below. 
With this amendment, Schwinger's result is 



G'miM) = Gi + Gn + Gni, (4.117) 

G, = ^!(P^, (4.118) 

E-T ' 

Gu = -A— ^7 ^ — , (4.119) 

2n^E-T{p- p'y E-T' ^ ' 



G 



III 



'2n^E-T 



dp / 1 
7 V4e?p(p-p')' + C(l-p)2 
1 



C 



E~T' 



(4.120) 



where 



r = P!,r'=P:!,e, = ^, (4.121) 

m m 2 

G = (e? + p^)(e? + p'^) , (4.122) 

M'i - 4m^ e? 
E = — = — i . 4.123 

4m m 

The prime on G'j^2 denotes the fact that we have subtracted the u = (1, 0, 0) part 
of G_x42 as required by the usual rules of second-order BSPT. Note that this Green's 
function is symmetric under p ^ p' and also {px,p'x) ^ {PyjP'y): symmetries that will 
be used in later simplifications of the integrand of SM^. 
is now 

2 

= Jd'pd''kd''p'd'k'{(Pi,0fi\k)V'"\k,S3,S,-p,Se,Se) 

Se^S— 

fGl + GlJ + GjJj\ (0). , , N/l/l, \ IA^C^A\ 

X W (p , Se, St, k , si, S2)(k |0i,o,o) (4.124) 



4m 

ml (/) + ml {II) + 6 Ml (III) respectively . (4.125) 
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Now we rewrite this in terms of hyperspherical harmonics and perform the integrations 



analytically. The variables are defined as 



[u = {uo,u)] ^ [ei,p] , [u = {u^,u')] ^ [ei,p'] 
[v = v)] ^ [ei,k] , [v' = «,v')] ^ [ei,k'] . 



(4.126) 
(4.127) 



See Appendix for a summary of the mathematical relations we use. The symbols 
appearing in Eqs. ([4.126|) and (|4.127| ) are explained there. Note that we use ei in 
these variable definitions, a choice that is completely general and turns out to be 
useful because we are taking expectation values of n = 1 states (e„ = ^ recall). The 
relations we use are 

4ef 1 



(k'|0i,o,o) 
1 



(ef + k'2)2 
(1 + ^;)^ 



Given these, 5M^ becomes 



(ef + k'2)2 Aef 
'el + k''f 



8e? 



el 



(1 + 



(4.128) 
(4.129) 
(4.130) 



ma 



327r2 



dVtpdVtpi dVt]~dVt}^i 



{E~T){Gi + Gu + Gni] 



(4.131) 



where 



S = ^t;(°)(k, S3,S4;p, Se,Se)i;^°np')'Se,Se;k',si,S2) . (4.132) 

Recall Eq. ([4.86|) for the definition of v^^\ Using the symmetries of the integrand, the 
sum over spins Se and can be performed and a simplification is seen to arise. The 
spin completely factors out of the momenta integrations. In other words, we have 

(4.133) 
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S = -(3^i + ^72)(p-p' + k-k'-2p-k') 



where 



gi = S1S3 + S2S4 , (4.134) 
g2 = 1 + S1S2 - S2S3 - 5154 + S3S4 + S1S2S3S4 • (4.135) 

Recall that Si = ±1, {i = 1, 2, 3, 4); i.e., the '|' has been factored out of these spins.Q 
So, in other words, instead of having to do sixteen twelve dimensional integrals be- 
cause the spin and momenta are coupled together, we just have to do one twelve 
dimensional integral that is independent of spin and then diagonalize the result in 
the 4x4 dimensional spin space with the spin dependence given by Eq. ( [4.133|) . 
We define the following integral 

where 

^ = J (l + n„)(l+<)2 [iE-T)iG, + Gn + Gni] 

x{p-p' + k^-2p-k')Y.^Y,{np)Y,inp,)Y:{nk)Y:,{nk>) , (4.137) 



and 

2 4 

8Ml = -^{3g, + g2)x. (4.138) 

For the quantities ^, x ^"^^ ^-^li labels /, // and /// imply the respective terms 
with G/, Gn and Gni above [see Eq. ( [4.117| )]. Also, the terms a, h and c above 
correspond to the respective superscripts in what follows. This integration will now 

be performed analytically. 

''^ In order to get these simple forms for gi and g2 and to see this spin/momentum decouphng it 
was useful to note the following simple relation: 5ss' ~ hs{s + s') [true because — 1]. 
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First the three Gi pieces. The mathematical relations used here are 

5\p-p') = ^±-^^s{%-Q,,)^^^^^5{n,-n,,), (4.139) 

p' = -^{1-Uo), (4.140) 
k-k' — > 3k,k'^, (4.141) 
p-k' ^ 3p,k'^. (4.142) 

Note that these last two relations are possible due to the rotational symmetry of the 
integrand. Then we expand these z-components of momenta upon the hyperspherical 
harmonic basis using the following simple relation (e.g. the Pz case) 

Now we recall the form of the hyperspherical harmonics (see the appendix on hyper- 
spherical harmonics for details), and their orthonormality and phase relationships 





= Yn,l,m{^) = fn,l{'^)Yl,m{OA) , 


(4.144) 


Yn,l,m 


— {~^y~^"^Y*^i_^ , Yi^m — (~l)™^/*-m , fn,l = 


(-1)'/;^ ,(4.145) 




= dfl = dQ^^^ dou sin^ u , 


(4.146) 


dQY^Y,, 


= St,^' , / du; sin^ u;fn,ifn',i = ^nn' , 





J dd"'Y*^Yv,m' = Sw6mm' ■ (4.147) 

After straight-forward application of these relations we obtain 

Att r , (1 — coso;) Btt^ 

• = — / dujsin^uj) { = 4.148 

Ci Jo (1 -I- COS a;) Ci 

• e = --£-fr^/vM/MM)'. (4.150) 
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For the Gu terms, we use the following relations 



1 1 m 



E-T' _£i _ 2e' 



^ - TT^^^ + O ' (4-151) 



1 

m m 



_J_ _ i2±^^2^„.,,^, , ,,,, 

These give 

{E-T)Gu = {l + u'^f{l + Uo)j:-Y,{n,)Y:{n,,). (4.153) 

^ n 

We use the rotational symmetry of the integrand and expand the integrand on the 
hyperspherical harmonic basis as was done for the three Gi terms. Then we have 

Svr^ ^ 1 / r duj sin^ u ^ ^ 



• 0"/ = 7^ E - / TT /2,i(^)/n,i(^) (4.154) 

2ei ^ n yjo l + coscj J 

• in = ^ E (/„ j (4.165) 

• a = -- E ^ f /' . (4.166) 

ei ^ n"' \Jo 1 + cos u J 

For the three Gjjj terms we use the same relations used for the three Gjj terms, 
and we use the rotational symmetry of the integrand to rewrite the appropriate pieces 
of the integrand in terms of 5^2, i,o a-s we did for the Gj and Gjj terms. However, we 
need to discuss one additional relation that allows the remaining 6M^{III) calculation 



to be done analytically. In Schwinger's 1964 paper [40] he gives the following formula 



11 OO 1 

T^n ^ = Ep""'-E^n,',™(^)^nV(^^'), (4.157) 
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where u and u' are of unit length and < p < 1.^^ Inside the brackets in Gm we 
have 



dp 



1 



p V4efp(p-p')2 + C(l-p)2 

1 dp 1 / 1 



(4.158) 



k p C \ {l - pf + p{u - u'Y 
Recall C = {e\ + p^)(e^ + p'^). Also recall that we are using the coordinate change 
of Eqs. (HH) and ( ^12^ ). Eq. (|Dl9| ) with e„ = d then applies and was used. In 
Eq. (I4.158D , < p < 1 and u and u' are of unit length, thus Schwinger's formula 
[Eq. (^4.1571) ] can be used and we have 

am{l + Uo){l + uy 



{E - T)Gin 



2ef 



fdp E 



n-2 



n 



1/^(1,0,0) 

Now, since n > 2 in this sum we can do the integral over p easily. 



dpp 



n-2 P 



n.-l 



n — 1 



1 



n — 1 



and we obtain 



{E-T)Ghi 



am{l + Uo){l + Mo 



/ \2 



1/^(1,0,0) 

For terms in ^ which contain Gni-, one obtains 



n[n 



—y,{%)y:{q. 



/// 



2 oo 

E 



(iu; sin^ u 



f2,l{^)fn,l{^) 



2ei n{n — 1) yJo 1 + costu" 
^■^This is easily derivable from a more general standard formula that Schwinger gives 



47r2 (w - u')^ ^ 2n 

^ ^ n— 1 '^-^ 



(4.159) 



(4.160) 



(4.161) 



(4.162) 
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37r2 ~ 1_ 



n=2 

37r2 °° 



/2,l(^)/n,l(^) 



/// 



E 



1 



1 + COS UJ 



n^[n — l) \Jo 1 + cos UJ 



f2,l{^)fn,l{^) 



(4.163) 
(4.164) 



ei n=2 

Now recall x = also notice that all the summands are the same, thus 

putting it all together we have 

1 



X 



3 3 ^ 

2 

n 



1 1 

n nin — V, 



1 1 



n^{n — 1) 



3 _ 3^ / I 



^^(n — 1)/ 
du; sin^ a; 



n 



(itusin u; 
1 + cos UJ 



f2Auj) fn,i{uj) 



/2,i(^)/n,i(^) (4.165) 



(4.166) 



10 l + COSCiJ 

The remaining sum can be done analytically. To see this, first define two integrals 



h = 



dVtp dVlpi 
1 + Mo 1 + m'o 



(47r) 



^ duj sin^ UJ ' 
1 + cos U) 



167r' 



(4.167) 
(4.168) 



1 + Mo 1 + < (p - p')2 3 
The last equality followed from rotational symmetry of the integrand. Thus I2 
327r'^/3. We can also calculate I2 a hard way which give^ 

r • 2 " 



167r^ 



4^ E - / TT /2,l(^)/n,l(^ 

^ n \ JO 1 + COS a; 



(4.169) 



Thus, we have 



00 -1 

E 



n=2 



n 



duj sin^ u; 
1 + cos UJ 



f2,l{^)fn,l{^) 



4 
3 



(4.170) 



Combining this result with Eq. ( ^.166| ) gives 



X 



3 /4 



8 V3 



We use Kr 



Px 



^2,i,i(^^p)], Eq. ( D.1S| ) and Eq. ( D.20| ), explicitly write the integrand in Eq. (4.168) out, and then use 
the phase and orthonormality relations of the hyperspherical harmonics to simplify the remaining 
expression. 



2 



(4.171) 

f [y2,i,_i(f7p)- 
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Thus, recalling Eq. ( [4.138| ), we have 



2 4 

5Ml = -^(3^1 + ^2), (4.172) 

where gi and g2 are given in Eqs. ( [4.134| ) and ( [4.135| ) respectively. This is the promised 
result written in Eq. (KTUDf). 



4.3 Singlet-triplet splitting: A simpler method 

Perhaps the most straightforward approach to calculate the singlet-triplet splitting 
is to just get busy and calculate, since the nonrelativistic Coulomb spectrum and 
states are so well known. This is exactly what is done in the previous section; however, 
as seen by the complexity of Subsection [4.2.3|, the calculation is complicated and at 



the level of a "Lamb shift calculation." Here we present a simpler method to calculate 
this shift. This simpler method uses a unitary transformation to "remove" V^'"' much 
in the spirit of Schwinger's early QED calculations |^ . 



First, we set up the notation. The exact eigenvalue equation is 

(ffo + V^)|$7v) =M^|$jv) , (4.173) 
where Mn is the mass of the state and 

• ($jv|$7V') = 5nn' (4.174) 

• 1 = X! J |P'SlS2)(pSiS2| = J |XS1S2)(XS1S2| 



Y,\^n){^n\ (4.175) 



N 



• (P'S3S4|V^|PS1S2) = V^(p's3S4;psiS2) (4.176) 

^^The idea behind this simpler method originated with Brisudova and Perry . To the consistent 
order in momenta required for fine structure, their transformation is equivalent to a Melosh rotation 
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(p's3S4|i?o|pSlS2) = 4(m + p )S {p - p')Ss,sJs2Si 



M^ = (2m + 5^)2. (4.178) 



m is the electron mass, —-Bat is the binding energy, and labels all the quantum 
numbers of the state. For notational purposes note that we label the final relative 
three-momentum with a prime, and that the initial and final electrons are labeled by 
"1" and "3" respectively, and the initial and final positrons are labeled by "2" and 
"4" respectively. V is given by Eq. ( |4.74|) , and its leading order and next to leading 



order matrix elements are given in Eqs. ( [4.85| ) and ( [4.89|) respectively. In zeroth order 
V is neglected and Eq. ([4.173|) becomes 



Ho\<Pn) = MI\<I)n) = (4m2 + AttiBn) \<Pn) ■ (4.179) 

This last equality defines our zeroth order binding energy, —Bn- Projecting this 
eigenvalue equation into momentum space gives 

(-Bn + —] Mp'sss^) = a / 7-^-^Mps3S4) , (4.180) 



m J 27r2 (p — p')^ 

the familiar nonrelativistic Schrodinger equation of positronium. 

After the simplification discussed below, as in the Coulomb gauge equal-time 
calculation, to obtain the ground state singlet-triplet splitting through order a^, only 
the wave function at the origin is required, which we thus record 



n is the principal quantum number, and / is the orbital angular momentum quantum 
number. This sets up our notation and we proceed with the simpler method. 
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The simpler method begins by acting on the Hamiltonian with a general unitary 
transformation with hermitian generator Q: 



JO) , ^(1) , ^(2) 



H = Ho + V ' + V ' +V" + 



H + t[Q,H] + -[Q,[Q,H]] + 



(4.182) 



(4.183) 



Now define Q by requiring its commutator with Hq to cancel V 



(0) 



v'°'+t[Q,Ho] ^ 0. 



(4.184) 



Putting this into Eq. ( [4.183| ) gives 



+ {V''' + V'^' + ■ ■ ■) e- 



+ 



iQ, ^iQ, V 
1 1 



(0) 



3! 4! 



iQ, iQ, tQ,V^°' +■■■} . (4.185) 



Note that H and H' have equivalent lowest order spectrums given by Hq (this can 
be seen easily by looking at matrix elements of the equations in Coulomb states, that 
is in states of Hq). However, the corrections to Hq in H start at order a^, whereas 
they start at order in H'. To summarize, we must diagonalize the following matrix 
element in spin space to consistently obtain the ground state singlet-triplet splitting 
in positronium through order a^: 



>1,0,0,S3,S4 



V" + 



(4.186) 



where Q is a solution to Eq. ( [4.184| ). Note that this is a first order bound-state pertur- 
bation theory shift. The quantum numbers are = {n,l,mi, Se, Se) — > (1, 0, 0, Sg, Sg) 
for the ground state. 
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In what follows we will solve Eq. ( ^4.184 ) for Q in the free basis in momentum 
spacej^ and then calculate the shift defined by Eq. ( [4.186| ). 



From the form of and Hq we see that Q has the following general form 



(p's3S4|zQ|pSiS2) = (5^(p-p')(p'-53-54K^|pSlS2) 



(4.187) 



where from Eq. (|4.184|) , R satisfies 
2m 



pS3S4\iR\pSiS2) — (p's3S4|ii?|p'siS2) . (4.188) 



Recall Eq. ( ^.86| ) for the form of v^^\ Thus, R is given by 



{pssSi\iR\psiS2) = ^'''^^fj'* {siPy - ipx) + ^'''J"^'^'* {s^py + ip^ 



2m 



2m 



. (4.189) 



Since Q is diagonal in momentum space it is a simple matter to calculate the 
contributions from Eq. ( [4.186|) . Define 



6m: 



0,0,S3,S4l^ |01,O,O,si,S2) 



■^l,0,0,S3,S4l 



iQ,V 



(0) 



,0,81,82) 



(4.190) 
(4.191) 



First, 6M' 
' 1 



2. 



d^p£p' {(pioo I p') (p 1 0100 ) (p'ss S4; psi S2 



(4.192) 



Using the rotational symmetry of the integrand, we can replace 



{P±-P'±y 



{Px - P'x? + {Py - P'yf + {Pz - P',f 



2 
3 



(4.193) 



(P - P')^ (P - P 

After this, the remaining integrals are trivial [recall Eq. (|4.181|) ] and we have 
bMl 1 11 

2^2q,4 ~ s\si^ s\sj> szsii, ~ ~^^s\'s-i^s-i's^ + s-£si,^ s\s-i^ s\sz ■ (4.194) 

^^This is the simplification: to solve for Q in the free basis; if Q is solved for in the Coulomb basis 
the (5M„^ calculation follows the one carried out in Subsection 4.2.3. 
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Next, 5M'^: 

' 2 



5M^ 



^l,0,0,S3,S4| 



iQy 



(0) 



^l,0,0,si,S2/ 



(4.195) 



1 f (0) 

'^^P'^V(01OO|P')(P|01OO) ((p's3S4|i-R|p'SeSe)(p'SeSe|V" ° |pSiS2) 



(p's3S4| IpSeSe) (pSeSe|i-R|pSiS2) ) • 



(4.196) 



Recalling Eq. ( [4. 85] ) and using Eq. ( |4.188| ) we have 



5Ml = ^1 rfW(0ioo|p)(p|0ioo) (p^p.p 



where 



F = ^(pSeSe|ii?|pSiS2)(p's3S4|f(°)|pSeSe) 



(4.197) 



(4.198) 



= 9 i{pSeSe\iR\pSlS2) - {p' SeSe\iR\p' S1S2)) (p's3S4 |f | pSeSg) , (4.199) 

using the fact that v^^^ is odd under p < — > p' in this last step. Using Eq. ( |4.188| ) this 
becomes 



1 
4m 



W^°^(p's3S4; pSeSe)v'^^\p'SeSe; PS1S2) . 



(4.200) 



Using the even symmetry of the rest of the integrand under the operations [p^ — * 
~Px,p'x — ^ ~P'x) {px < — > Py,p'x ^ — ^ Py) thls sum cau be simplified with result 



24m 



(3^1 + 92) (p - P 



f\2 



(4.201) 



where 



gi = S1S3 + S2S4 , (4.202) 
5f2 = 1 + S1S2 - S2S3 - S1S4 + S3S4 + S1S2S3S4 . (4.203) 
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Recall that Sj = ±1, {i = 1,2,3,4); the '|' has been factored out of these spins.0 



The result was written in this form to show the equivalence with Subsection [4.2.3 
Note how much simpler the SM^ calculation is in the present section. Putting it all 
together we have 

SM^ = --£^(3g^ + g2) I d'pd'p' {<P^oo\p'){p\<Ploo) (4.204) 
(3(71 + (72), (4.205) 



2 4 

ma 



24 



using Eq. ( |4.181| ) in this last step. 



Combining the results we have 

6M^ + 6M^ 1 , , 1 . . 1 



-^(3(71 + (72), (4.206) 



which is the same as Eq. ( [4.101| ) as was to be shown. 



In order to get these simple forms for gi and g2 it was useful to note the following simple 
relation: Sss' = ^s{s + s') [true because = 1]. 
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CHAPTER 5 



LAMB SHIFT OF HYDROGEN 



Experimentally the Lamb shift was discovered by Lamb and Retherford in 1947 
4^ . Later that year, Bethe submitted his seminal theoretical paper |Q. And quantum 



field theory began to look complete. For a review of some of these early calculations 



see Sections 19 and 21 of Bethe and Salpeter's classic text [38| and references within. 



Some selected early papers can be found in The agreed upon result — that has 
stood the test of time — that arose from this very active period, including QED effects 
through one loop, for the n = 2 and j = 1/2 levels of hydrogen is p6| 



5E. 



2Pi 



a^Ryd 

a^Ryd 
37r 



In 



m 



./5(2,0), 



ln2 + 



1 



6 8 5 



In = 



Ryd 



./5(2,1) 



(5.1) 
(5.2) 



where the last terms (+| and — | respectively) are the anomalous magnetic moment of 
the electron contributions and the — is the vacuum polarization contribution. Each 
of these terms has an interesting history The notation (3{n,l) will be explained 
later in this chapter; it is a particular average excitation energy of hydrogen that 



comes up in the calculation. Putting in the experimental parameters |^8| gives the 
following for the theoretical Lamb shift of hydrogen through one loop in covariantPI 



Actually, this "covariant" restriction can be removed as KroU and Lamb 
Weisskopf M showed. 



and French and 
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QED equal-time calculations 

^25i -^2Pi = 1052.19 MHz (27r;z) . (5.3) 

As motivated in the Introduction, we will only calculate the dominant part of this 
shift. Two modern experimental results are |^ 

^2Si - ^2Pi = 1057.845 ± 0.009 MHz (27r^) (5.4) 



and |2| 



^25i - ^2Pi = 1057.851 ± 0.002 MHz {luTi) . (5.5) 
For a modern status report on the theory of the hydrogen Lamb shift see and 



references within. For a modern review and text see 54 and references within. For 



some further selected references on the subject over the years see Here we must 



apologize for the inadequate references to the many papers on this subject. 

We proceed with an overview of our Lamb shift calculation. In hydrogen there is 
a small amplitude for a bound electron to emit and re-absorb a photon, which leads 
to a small shift in the binding-energy. This is the dominant source of the Lamb shift, 
and the only part of this shift we compute in this chapter. This requires electron 
self-energy renormalization, but removal of all the bare cutoff A dependence requires 
a complete 4th order calculation, which is beyond the scope of this dissertation. We 
work with a finite bare cutoff: A = m^j2. However, we do show that our results are 
independent of the effective cutoff, A,0 which validates our adjusting the effective 
cutoff into the range ma^ <^ A <^ ma, which is necessary to obtain the results in the 
few-body sector alone. 

= A — m — mp and A = A — to — m^: convenient definitions for discussing restrictions on 
typical binding-energy scales. 
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The energy scale for the electron binding energy is ma^, while the scale for photons 
that couple to the bound-states is ma. This energy gap makes the theory amenable 
to the use of effective Hamiltonian techniques. For simplicity, we use a Bloch trans- 
formation to remove the high energy scale {i.e., ma) from the states, and an 
effective Hamiltonian is derived which acts in the low-energy space alone. This ef- 
fective Hamiltonian is treated in BSPT, as outlined in Chapter ^. The difference 
between the 2S'i and the 2Pi energy levels, which are degenerate to lowest order, is 
calculated. 

We divide the calculation into two parts, low- and high-energy intermediate photon 
contributions. The low-energy photons satisfy |k| < A, while the high-energy photons 
satisfy A < |k| < m. A is the effective cutoff for the theory, which is chosen to lie in 
the range ma^ ^ A ^ ma. This choice lies between the two dominant energy scales 
in the problem and allows us to avoid near degeneracy problems. When an actual 
number is required we use 

A = aAj/a m ~ 6 X 10"'' m . (5.6) 

Note that the spectrum of the exact effective Hamiltonian is independent of A, but 
our approximations introduce A-dependence. The range for A is chosen so that this 
independence can be derived consistently order by order in the few-body sector alone. 

One further introductory comment, the high photon energy (A < |k| < m) part 
of the shift is further divided into two regions, A < |k| < 6 and 6 < |k| < m, where 
b is an arbitrary parameter chosen in the range ma <^ h m. This simplifies the 
calculation with appropriate approximations being used in the respective regions. 
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The result must obviously be independent of this arbitrary division point b, and is, 
unless "non-matching" approximations are used in the respective regions. 



We now outline this chapter. In Section |5.1| we review the basic theoretical frame- 
work of this light-front Hamiltonian approach for this Lamb shift calculation, and 
then in Section |5.2| we proceed with a discussion on the origin of the Coulomb inter- 



action in the electron-proton system. Section contains the heart of the Lamb shift 



calculation. In the final section. Section ^.4| , we summarize and discuss our results. 



5.1 Theoretical framework 

The proton will be treated as a point particle. The Lagrangian for the electron, 
proton, and photon system is (e > 0) 

C = -^F^uF^" + ^ + e^-m)ij, + i^p{t ^-e4-mp)ijp . (5.7) 

The reduced mass of the system is defined in the standard way 



nir 



m(l-m/mp + C(l/mJ)) . (5.8) 



m + nip 

Note that we take the limit nip/ m — > oo because we are only interested in the dom- 
inant part of the Lamb shift, but we will keep the reduced mass with a finite proton 
mass for the derivation of the Coulomb potential. The Lagrangian leads to the canon- 



ical Hamiltonian H in the light-cone gauge = that was derived in Section ^ 



The final form of the canonical Hamiltonian is written in Eqs. ( |2.49|) - (p.57|) . For a 



summary of the light-front conventions used for hydrogen see Appendix 0. 

Given the canonical Hamiltonian H we cut off the theory by requiring the free 
energies of all states to satisfy 



2 . \2 



Si < —pT— ' (5-9) 
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where A is the bare cutoff, and V = (V~^, P"*") is the total momentum of the hydrogen 
state. Then, with a Bloch transformation we remove the states with free energies 
satisfying 

p+ < < p+ ' (5-10) 
where A is the effective cutoff. The result is an effective Hamiltonian, Hx, acting 



in the low-energy [si < ^ ^1^:^ ) space. We do not discuss the derivation of Hx any 



further, but instead refer the interested reader to Section |375. 



Given Hx-, we then make the following division 

Hx = no + {Hx -Ho) = Ho + V , (5.11) 

where Tio is an approximation that can be solved nonperturbatively (and analytically 
for this QED calculation) and V is treated in BSPT. For this Lamb shift calculation, 
we treat the Coulomb interaction between the electron and proton to all orders in all 
Fock sectors. The test of Tio is whether BSPT converges or not and closely related: 
Is the A-dependence of the spectrum weakened by higher orders of BSPT? 

5.2 Lowest order Schrodinger equation 

The primary assumption we make in this QED bound-state calculation is that 
the Coulomb interaction dominates all other physics. After this assumption, the 
kinematic length scale of our system is fixed. 



1 1 137 



p ma m 
which then fixes our dynamic time and length scale, 

1 1 1372 



t 



p^/{2m) ma'^ m 
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as is well known, dynamic changes occur very slowly in this system. Note that in this 
QED calculation we will treat photons as free since they carry no charge and interact 
very weakly at low energies. After choosing Tio, the a-scaling of our BSPT is fixed, 
and the spectrum is then calculated to some desired order in a and a ln(l/a). 

In the Coulomb gauge, the Coulomb interaction appears directly in the canonical 
Hamiltonian, which of course is not true in the light-cone gauge. In the light-cone 
gauge, the Coulomb interaction arises from the leading terms in a nonrelativistic ex- 
pansion of a combination of two types of interactions in our effective Hamiltonian: 
instantaneous photon exchange and the two time orderings of dynamical photon ex- 
change. This is illustrated in Figure 5.1. These interactions arise from first and 
second-order effective interactions respectively. See Eq. ( |3.73|) for the form of the 
effective Hamiltonian, Hx. Now we discuss the derivation of the Coulomb interaction 
in the light-cone gauge for the electron-proton system. 

The time-independent Schrodinger equation in light-front coordinates that the 
sum of the three time-ordered diagrams in Figure 5.1 satisfies isQ 

X Vc (pNi^KSeSp) . (5.12) 

is the mass squared eigenvalue of the state (pN, where 'W labels all the quantum 
numbers of this state. The tildes will be notationally convenient below. We have 
introduced the following Jacobi variables 

= (xV+,K + xP^) , (5.13) 
p'^ = (x'P+ , k' + x'P^) , (5.14) 



50t 



For a derivation of Eq. (5.12) from the Schrodinger equation in Fock space see Subsection |4. 2. l| 



where this was done for the equal-mass case. 
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H H 





I q X [ 1 + 0(p/m) ] 



Figure 5.1: The effective interactions that add to give the Coulomb potential. "H" 
implies that the photon energy is greater than A. "L" implies that the electron kinetic 
energy is less than A. We choose ma^ <^ A <^ ma; these "H" and "L" constraints 
can thus be removed to leading order, and we are left with — 4e^/q^, the Coulomb 
potential. 



where and p'^ are the initial and final electron three-momentum respectively, and 

Pe +Pp = p1 +p; = P = {V+, V^) (5.15) 

is the total momentum of the hydrogen state. Note that k is a two-vector. The norm 
is defined by 

/ (i^K/(27r)^ / dx/{An) 4>*N{xKSeSp)4>Ni{xKSeSp) — 5nn' ■ (5.16) 

CO*' 
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Vc is the sum of the interactions given by the three diagrams in Figure and will 
not be written in all its gory detail. The leading order term of Vc in a nonrelativistic 
expansion is defined as Vc and is written below. 

The nonrelativistic expansion is defined in the following way. A coordinate change 
which takes the range of longitudinal momentum fraction, x G [0, 1] to Kz ^ [—00, 00] 
is defined 



l^z + \ K + K,i + 

X = ^ ; (5.17) 



V 

This step can be taken for relativistic kinematics, but there may be no advantage. 
Note that this definition of corresponds to the equal-time z-momentum of an 
electron in the equal-time center-of-mass frame Pe + Pp = 0. This intuition will guide 
us later. Then, the nonrelativistic expansion is an expansion in |p|/m; i.e., we assume 

IpI < m , (5.18) 

where we have defined a new three-vector in terms of our transverse Jacobi variable 
K and our new longitudinal momentum variable n.^ (which replaces our longitudinal 
momentum fraction a;) 

P = («:,/«.). (5.19) 

Note that the free mass-squared in the Schrodinger equation, Eq. ( |5.12| ), after this 
coordinate change, becomes 



X \ — X 



H — = V"^ +P +v"^o + P 



51r 



The interested reader should consult Eq. ( 4.16 ) and the discussion leading up to Eq. ( 4.39 ). In 



Chapter ^ the discussion is for the equal mass case, but is readily generalized to the unequal mass 
case using the rules in Appendix Note that in Chapter || we use a G-W transformation instead 
of a Bloch transformation; the Bloch transformation is chosen for the current chapter because of its 
simplicity. 
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(m + nipY + 2(m + nip) 



2m.r 8 rrir m?'m'^ 



+0 



iPl 



(5.20) 



which is invariant under rotations in the space of vectors p. is the reduced mass 
given in Eq. (|5.8| ). The longitudinal momentum fraction transfer in the nonrelativistic 
expansion becomes 



X — X 



Pz-Pz 



+ 



p^ — p'^ I ( m — m„ 



m + nip 2mnip{m + mp) 



O 



iPl 



m?{m + nip) 



The leading-order term of in an expansion in |p|/m is contained in 



(m + mp) I -— + —o-r^ii I '^-es^'5s,s;,^L , 



where 



q 



P - P 



A' - V"^' + P' + \/^^ + P 



X e(^^ ~ (^V^m2 + p'2 + y^m2 + p'^ 



(m + mp)^ + 2(m + nip) 



2|g. 



(m + nipY + 2(m + nip) 



(5.21) 



(5.22) 



(5.23) 



(5.24) 



(5.25) 



Note that 9^ and Oh are the constraints that arise from the Bloch transformation. 

It is convenient to define new cutoffs which subtract off the total free constituent 
masses of the state 



A = A — (m + nip) , 
A = A — (m + m„) . 



(5.26) 
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In the limit rUp — > oo we require A and A to be held fixed. Note that this implies 

- (m + rUpY _ ~ A^ (mp-^oo) ~ 

— rTT ^ — A -|- -— r > A , 

2{m + mp) 2[m + Trip) 

+ = A+_iL_(-Z^-)A. (5.27) 

2(m + nip) 2{m + nip) 

In terms of these new cutoffs, 9l and Oh above become 

To see the Coulomb interaction arising from the \ep) sector alone, we make the 
following requirements (which are motivated from the previous two equations) 



IP 



2 ^ 

< A -C IpI and A > |p| , (5.30) 



2m, 

also demanded for |p'| of course. These constraints will be maintained consistently 
in this chapter. Given these restrictions we have 

« 1 , (5.31) 
Oh ^ I. (5.32) 



Vc becomes 



K ~ (5.33) 



where 



K = (m + mp)' \ + I 5s^s'^5s^s'^ 



-(m + nip)'' ( — j 5,,45,p4 . (5.34) 
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This is summarized in Figure 5.1 . 

To finish showing how the Coulomb interaction arises in a hght-front Hamiho- 
nian approach, we need to know the Jacobian of the coordinate transformation of 
Eq. ( CT) , 



1 + 



m nir 



(p^ + 2kI 

2mm„ 



O 



iPl 



. (5.35) 



It is also convenient to redefine the norm 



Snn' 



= ^ J d-^p (plf{pSeSp)(f)N'{pSeSp) . 

S^Sp 

In this last line the tildes are removed from the wave functions by defining 



(5.36) 



')N[pSeSp) 



'Ml 

167r3 



[pSeSp) . 



(5.37) 



Putting it all together, the leading order expression for Eq. ( |5.12D in an expansion 
in I p I /m given the restrictions of Eq. ( |5.30| ) is 



P 



V(p'sX) 



a 



^n[ps'X) , 



(5.38) 



^ 2m,; '"^P' 27r2 J (p - p')^ 

which we see is the nonrelativistic Schrodinger equation of hydrogen, is the 
reduced massQ and — /?„ is the leading-order binding-energy defined by 

Ml - (m + 



2(m + rrij- 



(5.39) 



^^In what follows we set irir to its infinite proton mass limit, m. We kept above to show that 
it arises in Eq. ( ^.38 ) as in the equal-time case. 
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where A^^ is the leading-order mass-squared.0 The well-known bound spectrum is 

A. = (5.40) 
where Ryd = ma'^/2 of course. Note that Eq. (|5.38|) fixes the a-scahng of |p| 

|p| ~ ma . (5-41) 



Thus we see that the restrictions of Eq. (|5.30|) become 



ma^ -C A -C ma and A ma , (5-42) 

which will be maintained consistently in this chapter. 
5.3 Lamb shift calculation 

Given our leading-order spectrum, we proceed with BSPT. As advertised, this 
will be divided into low and high intermediate photon energy calculations. Before 
proceeding with these respective calculations, we discuss whether Coulomb exchange 
can be treated perturbatively or nonperturbatively in the respective regions. 

For the low-energy intermediate photon, the Coulomb interaction between the 
intermediate electron and proton must be treated nonperturbatively, whereas this in- 
teraction can be treated perturbatively for the high-energy intermediate photon con- 
tribution. This is seen by noting that each additional Coulomb exchange contributes 
a Coulomb matrix element and an energy denominator which is dominated by the 
larger photon-energy scale. Thus each additional Coulomb exchange contributes 



^ < ^ > (5-43) 



^•^The full mass-squared M|,, and the full binding-energy, —B^, are related by = (m -f- nip + 
Bn?. 
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where we used the virial theorem {(j)i^\a/r\(j)i^) = 2 Qma^/n^^. For the low photon- 
energy contribution, in principle |k|min = 0, and each additional Coulomb exchange 
can contribute and therefore must be treated nonperturbatively. Of course, 

when the Coulomb interaction is treated nonperturbatively, low-energy intermediate 
protons and electrons form bound-states from which long-wavelength photons de- 
couple. This nonperturbative effect leads to |k|mm ~ 16.64 Ryd; see Eq. ( |5.119|) 



below. For the high photon-energy contribution, |k|mm = A and from Eq. (|5.6|) each 
additional Coulomb exchange thus contributes at most 



8.5 X 10-2 , (5.44) 



A 

and can therefore be treated perturbatively. 

5.3.1 Low photon-energy contribution 

The low-energy shift arises from two sources which are shown in Figure |5]^. The 
first term comes from the low-energy photon emission part of the effective Hamil- 
tonian, (a| / (fx^dx~TCee'y\b) , treated in second-order BSPT, where the intermediate 
electron-proton are bound by the Coulomb potential (scattering states must be in- 
cluded too of course). Recall Eq. (|2.53|) for the form of TCee-y- The second term is the 



result of renormalizing the one-loop electron self-energy: a counterterm is added to 
the second-order self-energy effective interaction in {a\Hx\b), which results in a finite 
(except for infrared divergences) shift to the electron self-energy. This is shown in 
Figure |5.3| . The counterterm is fixed by requiring the electron self-energy to evolve 
coherently with the cutoff. The details follow those in Subsection [4.1.1| and will not 
be repeated here. 
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diagram L1 



diagram L2 



Ryd 

~ 3 71 



In ^ 



X 



P(2J) 



Figure 5.2: Low photon-energy contribution to the Lamb shift. Diagram LI represents 
the shift arising from treating photon emission below the cutoff A in second-order 
BSPT, where the intermediate electron-proton are bound by the Coulomb potential 
(scattering states must be included too of course). Diagram L2 is an effective self- 
energy interaction plus counterterm (shown in Figure |5.3| ) treated in first-order BSPT. 
/3(2, /) is the average excitation energy of the n = 2 levels; see Eqs. ( |5.113| ) and ( p. 1141 ) 
and the discussion above them for details. 



Note that the term where the proton emits and subsequently absorbs a photon is 
down by two powers of the proton mass with respect to the term where the electron 
emits and absorbs a photon. This result is subtle though, because it is true only after 
the light-front infrared divergences have canceled between two diagrams analogous to 
the ones in Figure ^12] . In Hppj: the |^ ■ A-^ term leads to infrared divergences that 
cancel in the difference, the a-d^/d~^ terms are the mentioned terms down by a power 
of rrip, and the mp/d^ terms flip the proton's helicity and thus do not contribute to 
this self-energy shift. 
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Figure 5.3: The sum of an effective self-energy interaction (arising from tfie removal 
of photon emission above the cutoff A) and a counterterm. The counterterm is fixed 
by coupling coherence as in Subsection [4.1.1| . The result is the interaction in diagram 
L2 of Figure O. 



Before proceeding with the calculation, recall that the exact binding-energy of 
hydrogen, —B^, in terms of its mass-squared is 

= {m + mp + BnY ■ (5.45) 

Multiplying out the right-hand-side, note that this implies 



2 



+ y = + (6.46) 

2[m + mp) 2{m + rrip) 

Recalling Eq. ( |5.39| ) , which is the exact definition of the leading-order binding-energy, 
—Pn, in terms of the leading-order mass-squared, A^^; and also defining the mass- 
squared corrections, SMff, in the form in which they appear in our calculation 

^ Ml + SMl, (5.47) 
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combined with Eq. (^.46 ), gives 



2[m + mp) 2[m + nip) 

In what follows, we will show that as far as the fermion masses are concerned, ~ 
mrrip and B"^ ~ consistently. Taking this as given for now, as rUp — > oo the 
leading term of Eq. ( ^.48| ) is 

Defining the binding corrections, —6B]\f, by 

Bn = Pu + SBn, (5.50) 



combined with Eq. (^.49|) , gives to leading order in 1/ 



rrir. 



a useful formula to be used below. This formula is useful because 6M'^ is calculated 
below, but 6Bn is the quantity that is measured. 

The low-energy calculation proceeds as follows. The first term of Figure 



IS a 



second-order BSPT shift which contributes the following to the mass-squared eigen- 
value [recall Eq. ( p.56| )] 



e 



Ll 



where k and s-y are the photon's three-momentum and spin respectively, and V = 
V^,V^^ is the total momentum of the hydrogen state t/^at, with quantum numbers 
= {n,l,mi, Se, Sp) = (z/, Se,Sp). 7^667 is the photon emission interaction given in 
Eq. ( p.53| ), where we have defined Vee-y = J (fx^dx^Hee-y- 6li restricts the energies of 



the initial, intermediate and final states to be below the effective cutoff The 
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explicit form of these restrictions is discussed below. Continuing the description of 
Eq. (|32|), 

r _ r (Pk^dk+e{k+) _ r d^k 

Jk ~ I ~ J (27r)3(2|k|) ■ ^ ' 

The last step comes from recalling that for an on mass-shell photon in the forward 
light-cone k'^ = k^ + k^ = |k| + A;^. The denominator and volume factors are 

Vol = {iJN{V)\^N{V))/V^ = IGtv^ 6^{V-V) = J d^x^dx- , (5.54) 
The two-body states are 

y<MPePp)bl{Pe)Blip,m , (5.56) 
\^N'{V-k)) = f ^Jktk^lQ7T^6^ {V-k-ki- fca) 

J kik2 

x^,,{k,k2)bl,{k,)Bl,{k2)\0) , (5.57) 



where (pN are solutions to Eq. (|5.38|) , the nonrelativistic Schrodinger equation of 



hydrogen, and 0Ar is related to 0Ar by Eq. ( p.37| ). 
Straightforward algebra leads to 



SM 



2 

LI 



X 



[l6n'6%k + k3- pe)] [l6n^6^{k + ki - p'j] C {p'„ V - p. 



X 0,,(A;i,P-A;-A;i)C,(A;3,P-A;-A;3)0.(Pe,P-Pe)^^, (5.58) 

where recall that v and v' are shorthands for (n, /, m/) and (n', m^) respectively, the 
usual principal and angular momentum quantum numbers of nonrelativistic hydrogen 
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(the "c" on the sum emphasizes the fact that the continuum states must be included 
also). See Eq. ( |5.55| ) for DENi. Ni is given by 



which after some algebra becomes 



(5.59) 



(Ana) 



2m^ — - i 



Tl' + t + 



+ 



^ k+ kt 

where we have defined a new object, 

p>{s) -- 



2ie _ k\{s,) _ p';{s,) ^ (2ie _ pjis,) _ 

k+ 



Pe 



p> + i s ejk p^ . 



ki 



,(5.60) 



(5.61) 



Notation: j, k = 1,2 only, s = ±1 only, s = —s, eu = —€21 = 1 and en = €22 = 0. 

We now discuss 9li and then simphfy 6Ml^ further. Recall Eqs. (|5.2CI| ), ( ^.26|) 
and (|5.27| ). We see that after the coordinate change defined by Eq. (|5.17|) , in the 



rUp — > 00 limit. 



e,, = e{x-T)e{x-T')e(~x-^ + oia'))e(~x-^ + oia') 



where 



T = |k| + -i/(p — k)^ + — m 



(5.62) 



(5.63) 



and T' = T\p ,p/. We have used the fact that the wave functions restrict |p| ~ ma. 

Recall that we are always assuming ma^ <^ A ^ ma. Thus, 6li can be simplified 



OfA ^9rX-T]9rX-T'] . 



(5.64) 



From the form of Eq. (|5.63|) , we see that this constrains the photon momentum to 
satisfy 



|k| < A 



(5.65) 
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to leading order in a. 

Note that the constraints coming from ^^i, summarized by Eq. ( p.65| ), require the 
photon three-momenta in bM^^ of Eq. ( [5.58| ) to satisfy 



k < pe,Pe- (5.66) 
Thus, Eq. ( |5.58|) can be simphfied further 



X (t^u' {p'e^'P - P'e) K' i.Pe.V - Pe) <Pu {Pe,V - Pe) ^^J^^ , (5.67) 

with the constraints ks = Pe, ki = p'^ and |k| < A. 

In the rup — > oo hmit, — > rup, and DENi becomes 

DEN^ = 2mp{Pn-Pn' -\k\)[l + 0{l/mp)] , (5.68) 

where we recalled Ai"^ = (m + m^)^ + 2 (m + nip) f3n and have used + k^ = 2|k|, 
valid for an on mass-shell photon in the forward light-cone. — /?„ is the binding energy 
of nonrelativistic hydrogen with numerical value Ryd/n^ for the bound-states. 

In the region of integration |k| < A = ma^/a <^ |p|, so after the coordinate 
change of Eq. ( ^.17|) [recall Eq. ( ^.19| )] we have 



+ ^ — + 



Ak^ / _L ,_L p-^Pz p'^p' 



p±+p'±-^-^\+0 (a'V^) . (5.69) 

m. m l V / 



k^m \ m m ^ 

The rest of the integrand is even under k^ —k^, so these terms in the last line, 
odd in k^, do not contribute. 



Putting it all together, recalling Eq. ( |5.51| ), we have 
SMh a r (Pk 



P J/' I I 
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where we recalled Eq. (|5.37|) , the relation between 0Ar and (pN- This is infrared 



(A;+ — > 0) divergent, but we must add diagram L2 of Figure to get the total 
low-energy shift. 



As previously mentioned, Diagram L2 of Figure |5]^ arises from the sum of an 
effective second-order electron self-energy interaction and a counterterm defined such 
that the electron self-energy runs coherently. The result of this interaction is to add 
the following to the binding 

X (p) 0. (p) ^ ^ ^ . (5.71) 



■\/p^ + — Y (p — k)^ + — 
Given the constraint |k| < A ^ |p|, this becomes 

^Bl2 ^ / ^ ^ ~ I*"') / / (p') 0.' (pO 

X (p) 0. (p) ^ . (5.72) 

This is the famous subtraction that Bethe performed in 1947 ||^. In our approach it 
arose as a consequence of coupling coherence. 

5Bl2 is infrared divergent {k~^ — > 0) as is 6Bli. This divergence arises from the 
first two terms of A^^i (the ones independent of p and p'). Noting that 

1 / fc"^^ , , \ (fc+^O) fc"' 



k+ J 2k+ 

we have 

(fc+-+o) 



Pn - Pn' - |k| 
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2k+ 



and these infrared divergent contributions from the first two terms of A^i cancel, 
leaving an infrared finite shift, 



6Br = 6B 



Ll 



6B 



L2 



X C' (p) 4>u (p' 



(5.73) 
(5.74) 
(5.75) 



2\ a y-^c 
37 



9(\-\k\) d'p d'p'^:{p')<p,,{p' 



X (p) (p; 



P ■ pV 1 I 1 

m'^ I /?n |k| |k| 



(5.76) 



This last step followed after averaging over directions as dictated by rotational in- 
variance. 

Eq. (|5.76|) is easy to integrate, and our final result for the low-energy photon 
contribution is 



SBr 



2a 

2a 
3nm' 



■^'(/3„.-/3„)ln 



Pn' — Pn 

A 



\{^u\p\4>u')\^ 



Pn' — Pn 



\{<Pv\i»\<P. 



(5.77) 
(5.78) 



where in this last step we recalled A ^ ma^. Note the A-dependence in the result. 
This will cancel after we correctly add the contributions coming from high-energy 
intermediate photons, which now follows. 

5.3.2 High phot on- energy contribution 

The high-energy shift arises from three sources which are shown in Figure 5.4. 
These are first-order BSPT shifts due to third and fourth order effective interactions 
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+ 



diagram H1 




diagram H3 




diagram H2 



a- 



Ryd_ / ln-b_ + In iH A 
71 V ^ ^ / 



Figure 5.4: High photon-energy contribution to the Lamb shift. These are third 
and fourth order effective interactions treated in first-order BSPT. These effective 
interactions arise from the removal of photon emission above the cutoff A. '6' is an 
arbitrary scale, required to satisfy ma <^b <^'m, that was introduced to simplify the 
calculation. Note the 6-independence of the result. The total contribution is a sum 
of Figure |5.2| and Figure p.4|. Note the A-independence of the combined result. 



see Section |3]^). The net result of these three diagrams is 



a a 



27r2q2 27r2q2 



:i + SVh) , (5.79) 



^"^In practice we calculate the effect from the third order effective interaction with instantaneous 
photon exchange and then generalize — ^^"^^j 6Vh — > ~ 27r°q^ ^Vh- It is assumed that the leading 
order helicity conserving term in a nonrelativistic expansion of fourth order effective interactions 
with dynamical photon exchange combine with the instantaneous photon exchange resulting in the 
Coulomb interaction in the same manner as what led to the original Coulomb interaction — the result 
is obvious, but the algebra was not carried out. 
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where q is the exchanged momentum of the electron, and 



5Vh = 6Vhi + 6Vh2 + SVh3 , 



(5.80) 



with 



5Vhi = \j^e(p'^^ -k+)e(p+-k+)NHieHi 



X 



{V- 



{V- 



6V, 



H2 



{M^ - M2) (Mq^ - M' 
^ 1^ 9 {p'+ - k+) 9 [pt - k+) Nh2 9h2 



,(5.81) 



5Vi 



H3 



'I Ik ^ (^'^^ " ^'') ^ " ^^^^ 

^ 1 1 

(Mo^ - M2) (M^2 _ j^^2 



(5.82) 



(5.83) 



The factors ±| in front arise from the form of the Bloch transformation [see the third 



order effective interactions of Eq. ( 3.73 )]. The vertex factors are given by 



Nhi = (A''i)(^^^p^_^^^,g^p^_^) , 

Nm = (A^l)(A;i-»p,_fc,fc3^p,-fc,p',^pe) 



(5.84) 
(5.85) 
(5.86) 



where Ni was defined in Eq. (|5.6CI|) . The free state masses are given by 



Mn 



M 
M' 



p2 _|_ ^2 _|_ /p2 _|_ ^2 



P ' 



+ rn? + + ml , 



k| + Y (p - k)2 + + Jp2 + m2 



|k| + V(P' - k)^ + m? + + ml 



(5.87) 
(5.88) 
(5.89) 
(5.90) 
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The Bloch transformation constrains the free masses of the states. As discussed in 
the low photon-energy calculation, the "L" restrictions in Figure can be removed 



given A ^ ma"^. However, the "H" restrictions lead to important constraints given 
by the 9h factors above, which we now discuss. They constrain the free masses to 
satisfy [recall Eqs. 



A< M -m-mp<A, (5.91) 
A< M' -m-mp<A, (5.92) 



where M and M' are defined in Eqs. ( |5.89| ) and ( |5.90|) respectively. 



As already mentioned, for convenience of calculation, we will divide this high- 
energy contribution into two regions, A < |k| < 6 and 6 < |k| < m {region one and 
region two respectively), with ma ^ 6 ~ ra^fa <^ m. Recall, ma^ ^ A ~ may/a <^ 
ma. We now show how this division into these two regions arises as a result of the 
constraints of Eqs. (|5.91| ) and ( ^.92|) . 



In this first region, |k| <^ m, and Eqs. (|5.91| ) and (|5.92|) become 



A < |k| + -^^^ — |k| < b, (5.93) 



A < |k| + ^^^-— |k| < 6, (5.94) 



2m 
'-1 

2m 

which is as we have already stated (recall that we always assume oo and drop 

the 1/mp corrections since we are just after the dominant shift). 

The analysis of the second region is slightly more complicated because |k| S> ma, 
and near the upper limit |k| ~ m. Since |k| ^ ma in this region, Eqs. (|5.91|) and 



( ^.92| ) both become 



b < |k| + Vk2 + m2 - m < A. (5.95) 
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This is just a linear constraint, 

. / 2m + 6 



\2m + 26 

which, since we choose b <^ m, becomes 

6 < Ikl < 



A / A + 2m 

< k < - ^ 

2 \ A+m , 



(5.96) 



A /A + 2m' 



(5.97) 



2 V A + m / 

We do not deal with removing the initial cutoff A dependence. A full analysis of 
this dependence requires a complete 4th order calculation, which is beyond the scope 
of this dissertation. We cut off the photon momentum at the electron mass, and 
proceed. Note that from Eq. ( ^.97] ), this choice corresponds to A^ = 2m^. The point 
of calculating these high photon-energy contributions is to show that our results are 
independent of the effective cutoff A by a consistent set of approximations [which is 
one step beyond Bethe's original calculation but still not at the level of calculation 
resulting in Eqs. (|5.1| ) and (|5l^)]. 

Taking a sample denominator, in the first region 



(M^ - M^) = (Mo + M)(Mo - M) 

/P' „, ( 
!^ 2mr, k — 

^ I 2m ' ' 



2m 



2mp\k\ , (5.9^ 



and in the second region 



(M^ - M^) = (Mo + M)(Mo - M) 



2mp ^m 



Vk^T 



m" 



-2m, |k| + — 



(5.99) 



Using these previous formulas, including 
, after summing, become 



m-p as mp — > oo, Eqs. ( [5.81| ) 



sv' 



a q 



^2 ,1 



H 



An J-i 



dcos6 



d\k\ 



\ m . 



(5.100) 
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in the first region (the "prime" on SV^j signifies the first region), and 



^ An m? Jo 2tt J-i Jb k 



Ikl , 

1 + c„ — (1 + cos( 



m 



|k| p-*- X p'^ Ikl / Ikr 

-Cd 2 2Se "2 V \ J 

m g-L m 



(5.101) 



in the second region (the "double prime" on signifies the second region). In 
the second region since the photon momentum is not necessarily smaller than the 
electron mass, we have kept two terms in the |k|/m expansion of the integrand. 
In the (9(|k|/m) terms we have introduced two constants, c„ and q, which denote 
numerator and energy denominator corrections respectively. Hereafter we set Cn = 1 
and Q = 1, as given by the theory. 

This term in the last line of Eq. ( |5.101|) that is dependent on the electron helicity 
Se/2 is written for completeness, but it does not contribute to the Lamb shift as will 
now be explained. An expectation value of x in the 2Si state vanishes, and 

2 

in the 2Pi state vanishes after an average over m/ is taken as dictated by rotational 
invariance. Also dictated by rotational invariance is an average over = Se/2, and 
then the term in question proportional to nig obviously vanishes. For details on why 
rotational invariance dictates these averages to be taken see the end of Appendix 0. 

Note that in combining SVhi, SVh2 and SVhs, many cancelations occur; most 
noteworthy, each contribution is individually infrared divergent {k~^ 0), but in the 
sum the divergences cancel. These final equations are easily integrated, and we have 



= -^^Inli] , (5.102) 



2-7r \X 

a q fm 

T 1^ T 

2n \ b 



^y'u = -T^^^lnf^) . (5.103) 
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In the second region note that the (9(|k|/m) terms coming from numerator and energy 
denominator corrections cancel, leaving the 0{1) piece alone. The combined high- 
energy contribution is 



,V„ = SVl,^Sv;i^-^'^U.Cf), (5.104) 

which is independent of b, as required for consistency. Recall that q = p' — p: the 
difference between the final and initial electron momenta. 

From the definition of 6Vh [see Eq. (|5.79|) ], this correction shifts the energy levels 



an amount 

SBh = / d'p Sp'ct^l (pO ( (^^^ ) 0. (P) • (5.105) 

Combining this with Eq. ( |5.104[ ) gives 

where in this last step we averaged over directions and noted that the wave function at 
the origin is real. For more details on this averaging over directions see Appendix ^ 

5.3.3 Total contribution 

In this subsection we combine the results of the last two subsections for the low 
and high photon-energy contributions, and perform the required sums/integrations 
to calculate the total shift between the 25*1 and 2Pi energy levels of hydrogen. 
Adding Eqs. (|5.78| ) and (|5.107|) gives for the total shift 

A 



6B = 6BL + 6BH = ^-^Y.'{(3n'-(3n)\n 



\V\<\>V')'\ 
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For the second term we have 

2 



d p (j)u {p] 



{(2.)^J. = 0)y=^^(^p„. (5.109) 



The {2n)^ factor arose because of our normahzation choice [recall Eq. ( |5.36|) ]. 

The first term of Eq. ( |5.108|) is the famous Bethe log and must be calculated 
numerically, summing over all bound and continuum states. Following standard con- 
vention we define an average excitation energy, /3(n, Z), 



In 



Ryd 



-) if^n' - Pn)\{(j)u\p\(f)u' 

/ v' 



|p|0,,)l' • (5.110) 



The sum on the left is evaluated by standard techniques [He = /{2m) — a/r\, 

- Pn) K0.|P|0.')|' = \{<Pu |[p, ife] ■ P + P ■ [i^c, p]\<Pu) 



0, i-a/r) (f). 



p-,-2V(-a/r) 0, 



2a(^)%,o. 



(5.111) 



This vanishes for I ^ 0, but the average excitation energy, i3{n,l), is defined [it is 
just a way to catalogue the numerical sum on the right of Eq. ( [5. 1101 ), the quantity 
we need to know] with the sum on the left-hand-side set to its value for Z = 0. In 
summary, l3{n, I) for all states is defined by 



= E (/9n'-/3n)ln 



Ryd 



■ (5.112) 



Y Ryd I \ n 

Without further ado, this sum has been evaluated by R. W. Huff His results for 
the n = 2 levels are 



P{2,0) = 16.63934203(1) 
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(5.113) 



/?(2,1) = 0.9704293186(3) , (5.114) 

where the figures in parentheses give the number of units of estimated error in the 
last decimal place (R. W. Huff's estimates). 
Combining the results: 

= ^ln(=P-] , (5.115) 
r n 2a , / Ryd \ ^ ( ma \ ^ 

= ^JJyl_] . (5.116) 

Note the cancelation of the A-dependence. Thus, the dominant part of the Lamb shift 



IS 



SB, . ^ 6B.S. - 6B,P. - (5.117) 
= (1047 - 4) MHz {27ih) = 1043 MHz {27ih) , (5.118) 



where we use and the average excitation energies of Eqs. ( |5.113| ) and ( p.ll4| ) 



Note that the 2Pi shift is only about one-half of a percent of the 2Si shift (but both 
shifts are "up"). 

5.4 Lamb shift summary and discussion 

In a light-front Hamiltonian approach, we have shown how to do a consistent 
Lamb shift calculation for the n = 2, j = 1/2 levels of hydrogen over the photon 
energy scales 

ma^ A ^ ma ^ b ^ m , 
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Figure 5.5: The experimental n = 2 hydrogen spectrum: fine structure, Lamb shift 
and hyperfine structure. F = L + Sg + Sp. 



with the choices ma^ ^ A ^ ma and ma <^ b m. In a consistent set of diagrams 
we showed how A- and 6-dependence cancel leaving the dominant part of the Lamb 
shift, 1043 MHz. For completeness, the experimental n = 2 spectrum of hydrogen is 
shown in Figure 



If we compare the three regions we see the following results (we only compare for 
the 25*1 shift since the 2Pi shift is negligible within our errors) 



(0 < |k| < A) 5B^ 



a^Ryd 

371 



In 



A 



16.64 Ryd^ 



46 MHz {27rh) , (5.119) 



(A < |k| < b) 6B^^^^ r^^f^ln ~ 667 MHz {27rh) 



(b < Ikl < m) 6B, 



a^Ryd 



■In, , 

Sir \b 



m 



334 MHz (277/1) 



(5.120) 
(5.121) 
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where we used A = ma^fa and h — m^/a, consistent choices used throughout this 
chapter. As expected on physical grounds, the dominant photon momentum that 
couples to hydrogen is 

|k| ~ l/oo ~ ma . (5.122) 

That is, the photon wavelength that couples the strongest to the neutral hydrogen 
system is of order the size of the system. As seen above, the effects of photons of this 
momentum amount to about two-thirds of the Lamb shift, the dominant part of this 
experimentally observed shift. 
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CHAPTER 6 



SUMMARY AND DISCUSSION 



In positronium, in a consistent set of approximations, we analytically calculated 
the spin splitting for the n = 1 levels through order a^. To go beyond this, and 
calculate the spin splitting through order for example, it is clear that the effective 
Hamiltonian Hx must be derived through fourth order. 

Kaluza and Pirner were the first to numerically calculate the ground state spin 
splitting of positronium (neglecting annihilation channel contributions) through order 

in a light-front hamiltonian approach ,0 but they were forced to make ad hoc 
assumptions because their Hamiltonian depended on the full eigenvalue of the prob- 
lem. We avoided these assumptions in our approach, and performed the calculation 
analytically. 

The dominant part of the Lamb shift was calculated with result 1043 MHz. This 

turned out to be accurate, and within a consistent set of approximations was shown to 

be independent of our final cutoff A. It was important to establish this independence 

to validate the placement of our final cutoff in the range ma^ ^ A ^ ma — which is 

necessary to obtain the result in the few-body sector alone. We did not show that the 

result was independent of the initial cutoff A. We fixed A such that the photon energy 
A recent thesis that summarizes and extends this work is js^ . 
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in the respective loop satisfied |k|p/ioto„ < m, and did not derive the A-independence 
of the result. To obtain more precision, the full one-loop renormalization must be 
performed of course.^ A final note on the Lamb shift calculation: Each of our five 
diagrams (of Figures and were infrared {kpf^^^^ — > 0) divergent; however, 
both the sum of the two low photon-energy diagrams and the sum of the three high 
photon-energy diagrams were infrared finite. This is an encouraging result for future 
work. 

We close with a discussion of rotational invariance in the context of light-front field 
theory, and how the results of this dissertation fit into the big picture. The transverse 
rotation generators are dynamic and our regulator is not rotationally symmetric,^ 
thus non-covariant counterterms will in general be required to obtain a dynamic 
solution which adheres to the principles of rotational invariance. An initial study 
pO| at the one- and two-loop level in Yukawa theory showed in perturbation theory 
that rotationally invariant results follow with an appropriate renormalization of the 
kinetic and vertex masses. We continued the study of rotational invariance in this 
dissertation from the perspective of the QED bound-state problem. In Chapter § we 
showed in a consistent set of approximations through order in the ground state spin 
splitting of positronium that the result is rotationally symmetric and followed from the 
canonical structure of the interactions;^ however, this result could only be obtained 

after including the dynamic shifts from second order bound-state perturbation theory 

^^For an initial study of one-loop perturbative renormalization of QED in a light-front Hamiltonian 
approach with various regulators see ]59[ |. 

Since the generators are dynamic, we do not even know how to regulate in a rotationally sym- 
metric fashion until the dynamic solution is known. 

^^This necessarily included one loop renormalization of the fermion kinetic mass. Strictly speak- 
ing (from the transverse renormalization group sense), an infinite number of counterterms were 
required — these were fixed by coupling coherence. 
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which necessarily included a sum over all intermediate bound and scattering Coulomb 
states. 

As is no surprise, the generators are dynamic, and rotational invariance comes 
about in a dynamically complicated way even at the no-loop level for the bound-state 
problem — analytic understanding of this was presented in Chapter ^ and may be 
useful for future work. 

In Chapter ^ we obtained the dominant part of the Lamb shift of hydrogen con- 
sistently through order a^ln(l/Q;). This showed that the level classification — 251/2 
and 2Pi/2 for this example — as dictated by rotational invariance was maintained in 
this example. This calculation necessarily included one-loop electron kinetic mass 
renormalization which followed from coupling coherence. 



Neither example required the problematic diagrams mentioned in |Q with an 
instantaneous fermion line immediately preceding external self-energy shifts, thus 
we could set the vertex mass to the kinetic mass (actually it was required to get 
correct results). However, since the transverse rotation generators are dynamic, more 
nonperturbative analysis (for the positronium example) was required than in an equal- 
time calculation. Analytic understanding of the bound-state problem in QED in a 
light-front hamiltonian approach has been advanced. 
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APPENDIX A 



SUMMARY OF LIGHT-FRONT CONVENTIONS 



In this appendix we write our light-front conventions for the electron, proton, 
and photon system (including the antiparticles). The reader unaccustomed to light- 
front coordinates may want to read Chapter |] before consulting this appendix. Our 
conventions are 



ct = 
A+ 

V'e- 

A- 
A+ 



yO _|_ ^3 ^\^Qj.Q jg g^j^y four— vector. 



" 


' 




" -2i' 


2i 





; 7 = 






a' 
a' 



1 




z = 1, 2 ; cr' are SU{2) Pauli matrices 

; A- 




1 











e for electron, p for proton 







1 

L id+ 





-^{[a' {^^^-eA')+lmJ,]ip} 
—2 (9* 

2e(efe-elee) ; e>0, 
. 
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In momentum space the dynamical field operators are expanded as (at x'^ — 0) 



with 



Y: [ (e>,(g)e-«- + h.c.) , 

:=±1 

II Jv 



s=±l 



s=±l ''P 



Xi 



1 



1 ' 





1 



Note that s = —s. Also, we are using the shorthand 



-fip) 



I 2 , 2 



(27r)4 / / ^ v/-/ J i67rV 

The fermion helicity can only take on the values ±1/2, however we define hs — s/2; 
therefore, "s" can only take on the values ±1. The non- vanishing commutation (anti- 
commutation) relations and the free Fock states are given by 

[ax{q), aUq')] = 16nYS'{q - q')dxx> , [ d'{p) = S'{p^)5{p+) ] , 

{bs{p),bl,{p')} = {4(p),4(p')} = iQ7r-Ys'{p-p')Sss' , 

{BM, bUp')} = {D,{p), dUp')} = 167rV<^'(p - p')Sss' , 

{piSl\p2S2) = l6TT-^pt6^{pi -P2)5s^s2 : IPI^i) = &ii(Pl)|0) , etC . 
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APPENDIX B 



MATRIX ELEMENTS OF THE CANONICAL 
HAMILTONIAN USED IN THE DISSERTATION 



In this appendix, the matrix elements of the canonical Hamiltonian H in the basis 
of the free Hamiltonian h that were used to derive the results in this dissertation are 
given. See Eqs. ( p.49| )-( p.57| ) for the operator form of H. These matrix elements were 



worked out by Allen |T^. Note that these matrix elements correspond to normal- 
ordered interactions with zero modes dropped.]^ 

First, some comments and notation for clarity and completeness of this appendix: 
H = h + v where h is the free field theoretic Hamiltonian. h\i) = ei\i) and the matrix 
elements below are in this basis (e, p, e, 7) labels represent (electron, proton, 

positron, photon) respectively. Transverse momenta are written a.s p\, P2, ... {i = 
1, 2 only as always) and longitudinal momenta are written as pf, P2, .... Repeated 
transverse indices are to be summed. The 1/2 is factored out of the fermion helicities 
so all helicities below are written as si, S2, . . ., where Sj = ±1 only. Also note 
Si = —Si- m is the renormalized electron mass, rrip is the renormalized proton mass 

and — e [e] is the renormalized electric coupling of the electron [proton] (e > 0). The 

^^As explained in the dissertation, coupling coherence fixes all counterterms, and it is our con- 
jecture that it is not necessary to start with "self-inertias" or "zero-modes." If they are required 
for the renormalized Hamiltonian Hx to run coherently, they will arise through the process of 
renormalization. 
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photon polarization vectors are 



(B.l) 



We also use 



p^{s) = p'+isejkp'', 



(B.2) 



with ei2 = —621 = 1 and en = €22 = 0. A final convenient notation is 



a \v\ h) = 167T^d^{pa - Pb) {a\v\b) , 



(B.3) 



where Pa andp^ are total three- momenta of the respective state and S^{p) = d'^{p^)S{p'^ 
The eleven canonical matrix elements used in this dissertation are 



(e273|t^|ei) = -e^/: 



(es \v\ 6172) 



P1P2 



Pi Pi 
2i ( \ 

72 



Pi Pi 



eypipt 



(6273 1^^ I ei) = 



eypiPi 



V ^^3 Pi P2 , 

2i^_ (I l_ 

'-'^ 'A pt pt Pt . 



2i 



1 



-ex/PiPt 



.^""'^^ [%t Pt 
2i 



{P3\v\pi'y2) = eJpJpi 



pi pi(si) _ pi*is3) 
pt 

'''^'^Pt Pt Pt J 
^ (1 



(B.4) 



, (B.5) 



, (B.6) 



(B.7) 



(B.8) 
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(63 1^16172) = eJpJpi 



I ^ -L -L -±- 



Pi Pi 
2i 



+ ^sisa^sjs.—^Sim — r ir 



V2 



,P3 Pi 



(73 \v\ 6162) = -e\lptP2 



Pa Pi Pt ) 



2i 



1 1 



(6263 \v\ 71) = -esjptpt 



V2 \vX P2 

_ P^*(S3) _ P2is2y 
^Pl pt pt , 

2i f 1 1 



,P2 P3 



(e3P4 |v| eipa) = -4e^ 



2 V Pi ^2 



{pt - pt) 



2 ^51S3^S2S4 7 



(e3e4|w|eie2),^,;,„„^, 

{63.64 \v\ 6162) ^^^ihil 



(Pt-Pt) 



2 "S1S3"S2S4 ) 



(P]^+P2) 



2 "S1S2"S3S4 ) 



(B.9) 



(B.IO) 



(B.ll) 
(B.12) 

(B.13) 

(B.14) 



where these last two add as (6364 |^;| 6162) = (6364 It;] eie2),^^^„„g,+(e3e4 It;] eie2)„„„i^i;. 
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APPENDIX C 



4TH ORDER BLOCK AND 3RD ORDER G-W 



In this appendix we record the fourth order effective Bloch interactions as promised 
in Eq. (p.73|) of Chapter |^. Also we record the third order effective G-W interactions 



as promised at the end of Section |3]2 
First, fourth order Bloch: 



2 II(«I^AK)(^I^Ab')0'I^A|fc)(fc|^^A|&) 



1 



Yu{0'\'^A\i){i\vA\c){c\vA\j){j\vA\b) 



A-aiAajAak A^iAhjAhk 



1 



AaiAfij A(,j AfjiAfij A^i 
1 1 



^ai'—^ci'—^C'j 



\ II(«I^a|c) {c\vA\i) {i\vK\i) 0>a|&) 



{a\vA\i){i\vA\c){c\vA\d){d\vA\h) 
2 .^^^^ AuA^,Adi 

1 ^ {a\vA\c){c\vA\d){d\vA\i){i\vk\b) 

9 ^ AAA 

^ r r1 i '-^ai'-^ci'-^di 



AaiA^i 


Abj 


AaiA 


\ ' \ 




- ' \ 












--] 


A^Abj 1 


I A,, 


An J 



1 1 



AciAfij A(jj Acj 



(C.l) 
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And now we record third order G-W: 



_{3) 



E 



ab 



VAiaVAabVAbj 
+ 



6 



+ 



6 



+ 



jb'-^ia 



A, ft A 



+ U ^ 



J)' 



(C.2) 



where 



— f XabfAab^bai^abj^jab + f \aj f t^aj^ jai^ abj f Kab^baj t 
^2 = f XjbfAjbQjbiaQjbaQajb + f \aj f t^aj^ jai^ jbaf Kjb^bja j 
^3 = f XbjfAbjQjbi^iabjfAia^iab 1 
^4 = f Xbjff^bj^jbi^abjfAab^hai ■ 

This is for /aij = ^^(A— |Ajj|). Recall fx+fx = 1 Ajj = ^i— ^j. Note that the labels 
H, j, . . .' and 'a, 6, . . .' span the whole energy space [unlike the notation used for 
Bloch; see below Eq. ( p.73| ) for the notation in this case]. Qabcd = 6'(|Aafe| — |Acd|) and 



Qabc = 6{\Aab\ - lAfecl). Note that '(i < — > j)^' implies 'vAiaVAabVAbj 



VAibVAbaVAaj 



and (i < — > j) in the remaining pieces of the respective term'. As advertised we see 
that the denominators are constrained to be between A and A (A < A): there are no 
small energy denominators (as long as A is not too small). 
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APPENDIX D 



HYPERSPHERICAL HARMONICS AND FOCK 
COORDINATE CHANGE 



In this appendix we list some useful mathematical relations used in Chapter ^ 
The conventions followed here are the same as in Ref. [^. These hyperspherical 
harmonics are given by 



(D.l) 



where 



< \m\<l<n-l 



(D.2) 



These function labels v = {n, I, m) correspond with the standard nonrelativistic 
quantum numbers of hydrogen. These spherical harmonics, Yi^rn{(^,4'), are given 
by H 



Yi 



Lm \ 



2/ + l(/-m)! (-l)'" , ^.^ 



^ 47r (/ + m)! 2H\ ' dx^ 

where x = cos6. These other functions, are given by 



{x' - lye 



I inup 



(D.3) 



-I 



2nin-l-l)\ I 



\ 7r{n + l)\ 



(dcosu)' 



Cn-l(cOSUj) 



(D.4) 
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where C„_i(cosco') are Gegenbauer polynomials. The first few are |63 

Co{y) = 1, Ci(y) = 2y, C2(y)= C3(y) = 8i/='-%, 

C,{y) = IQy^ - I2y^ + 1 , C,{y) = 32y' - 32y' + Qy . (D.5) 

The phase and orthonormality relations of these above functions are 

Yn,l,m = , = (-1)'"F,:„„ , = (-1)7;, (D.6) 

and 

J dd^Wi*^Yi, = 6u'6mrn' , J duj Sm"^ UJ f*ifn' , I = ^nn' 

dnr^Y,, = 5,,, , (D.7) 

where dO. = rffi^*' = d^l^^' dcu sin^ u;, 0<u;<7r, 0<^^<7r and < (p < 27r. Combin- 
ing the spherical harmonics and respective derivative of the Gegenbauer polynomials, 
for the first few hyperspherical harmonics gives 

_ 1 v- _ v^coscj 

^1,0,0 - /— ^ 5 ^2,0,0 - 

—ie~'^^ sin u sin 9 —i \/2 sin cu cos 9 

^2,1,-1 = , ^2,1,0 = 

vr TT 

ie*"^ sin cij sin ^ /t» o^ 

^2,1,1 = • (D.o) 
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For further harmonics we refer the interested reader to Appendix 2 of Judd's text 
where this is done quite nicely. 

For the coordinate change in the Coulomb Schrodinger equation of positronium 
[see Eq. ( [4.39| ) and the discussion that follows it], for B^^ < 0, we define 

mB^ = -ej , (D.9) 
u = (mo,u) , (D.IO) 
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2 2 

uo ^ cosM = %^, (D.ll) 
e2 +p2 

u = 5111(0;) = 5111(0;) [5m(6') co5(0), sin{6) 5iii(0), co5 



^ 2enP 

Note that Mq + = 1. Conver5ely thi5 coordinate change give5 



(D.12) 



-u , (D.13) 



1 +M0 

el + v' = (D.14) 

We al50 have 

dap = sin^ ujdujdnf = (t^^^I d^P , (D.15) 



After thi5 coordinate change, the nonrelativi5tic Schrodinger equation of po5itron- 
ium i5 trivially integrable a5 de5crlbed below Eq. ( [4.39| ) in Chapter ^. Here we record 
the re5ult for the bound wavefunction5 and eigenvalue5: 



4(e 



(P) = TZ^Tl^YA^p) (D-17) 



and 



^ ma^ ma 



respectively. 

Finally, a most useful relation is 



ip-pf = M±p!M±p:!)|.-.f . (D.19) 

/|2 



This is useful because we can expand \u — u'\ as follows 



1 27r^ 

J2_Y,inp)Y:iap,). (D.20) 



\u — u'\^ ^ n 
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This completes the brief discussion on the hyperspherical harmonic mathematical re- 
lations used in Chapter ^; actually, as shown in Subsection ^.2.31 , the ^M"^ calculation 
requires some further formulas, which are given as they are needed. 
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APPENDIX E 



Ml -VS- IN POSITRONIUM AS AN ILLUSTRATION 



In this appendix we invert the equation 



Ml ^ {2m + B,)\ [El] 



and obtain the a-expansion for the binding energy, —B^. In Chapter ^ we set up a 
procedure to calculate the mass-squared of positronium. This gave 

Ml = M% + b^m^a^ + O [a^) , (E.2) 

with 

M% = Am^ + AniB^ . (E.3) 

For the leading-order spectrum of Hx we obtained 

^ 1 ma^ ,^ , 

= —A—T- E-4 

4 

Taking a square root of gives 

Recall that ^ = a^Ryd. Now, in this work, Eqs. (|4:T0^ )-( ^:T05| ) are the resuhs of 



our calculation of the spin splittings of in the ground state of positronium. These 

139 



were derived in the form of Eq. ( |bL).2|) with result 

7 

b^itriplet) — h/^{singlet) = - . (E.6) 

Given Eq. ( [E.5| ) this imphes 

,-B , , = -a^Ryd + O (a^) . (E.7) 

This we recognize as the well known result for the positronium system. A final note 
is that if the physical values of the fine structure constant and Rydberg energy (ysTo 
and 13.60 eV respectively) are applied to this previous formula, the result agrees with 
experiment to one-half of a percent [|65|1 . 
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APPENDIX F 



AVERAGING OVER DIRECTIONS 



We start this appendix by calculating the following Coulomb matrix element 



IP - p 



(F.l) 



to verify the step taken from Eq ( 5.106 ) to (|5.107 ). It is useful to define another 
integral 



IP -p 



'^2 



IP) 



(F.2) 



Now note that 



/ = lL + h = Jd'pJ d'p' C(p')0.(p) 



(27r)^ /ma 
TT \ n 



6, 



LO 



(F.3) 



where in this last step we recalled Eq. ( ^.109|) and the fact that the wave function at 
the origin is real. 

For / = 0, the wave function satisfies 



1,0,0(P) = 0n,O,o(|p|) • 



(F.4) 



Thus, by symmetry, for / = 0, 



1^ _ (2vr)^ /ma\3 
3 Sir \ n J 
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(F.5) 



and 



3 Sir \ n 

For / ^ 0, first note that / = 0. Thus, for / ^ 0, 

l± = -h ; (F.7) 

we will calculate Iz below which then implies I±. Next note that in position space 

/ = -27r2 j d'x 0: (x) (v'^^ 0. (x) , (F.8) 
using V^i^ = —AnS^ (x). Thus, for / 7^ 0, in position space 



X 



Note that there is no |x| ^0 ambiguity in this previous equation because for / 7^ 0, 
the wave function vanishes at the origin. Carrying out the derivative gives 



-27r^ / d'x (j^l (x) rrr^ W • (F-io) 



1X13 



This matrix element was performed in the first appendix of Bethe and Salpeter's 
textbook |6^. We use two of their formulas, (3.26) and (A.29).0 Eq. ( |F.1CI| ) integrated 



gives 



= -27r2 r-3 c(/, , (F.ll) 



with 



' \ (2/ + 3)(2/-i) ; ^ ' 

^"^Warning to the reader: In this text, they use atomic units, H — c — m — a — 1, som and a 
have to be placed back into the formulas. 
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Thus, recalling Eq. our result for Z 7^ is 

= 2n^V^ c{l, mi) . (F.14) 

For / = 1, J_L is not zero, so what is going on? The answer lies in the fact that we 
really want to take matrix elements in the |j, rrij) basis not the \mi,ms) basis,Q and 
based on rotational invariance our results are assumed to be independent of rrij. To 
proceed, note that we assumed is constant in our calculation of the dominant part 
of the Lamb shift since it is a self-energy shift. Also note that the result, Eq. ( [F.14| ), 



is even under mi — > —mi. Below we will use these facts to show that the Clebsch- 
Gordan coefficients (specifically we show it for the non-trivial case — the 2Pi states) 
imply 

{j = l/2,m,\V\j = l/2,mj) = 

11'^ 



Y^T^ {mims\V\mims) . (F.15) 



For now just take this result as given. Now note that I± given by Eq. ( p^'.14D averaged 
over mi vanishes, 

^i:/x=o. (F.16) 

mi=-l 

where we used 

1 ' o 1 



I E = -/(/ + !), (F.17) 



2/ + . , 

mi=-l 



an obvious result after the answer is known. This result [Eq. ( [F.16|) ] was used in the 



step that led from Eq. ( p.l06| ) to Eq. ( p. 107] ) in the dissertation. 



^^Recall throughout this dissertation we have specified electron hehcity by Se/2 = ±1/2. Here we 
set Se/2 — > rus- Also we use 's' to denote the electron spin, that is s = 1/2. 
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Now we show that given 'm^ = constant', 'the symmetry under rrii — > —ttii 
and 'independence of m^' Eq. ( [t'MSI ) follows. Using the standard Clebsch-Gordan 
coefficients [Q, the 2Pi expectation values are 

r 



1 1 
22 



V 



2 

1 1 

22 



1/1 
3 °2 



V 



1\ 2 / 1 
3 \ 2 



V 
V 



1 - 



1 

1 - - 
2 



fF.18) 



and 



1 1 

2 ~ 2 



V 



1 1 

2 ~ 2 



V 



- 



1 



2 

2^2 



2 
3 



-1- 



(F.19) 



where on the left of the equal sign the states are in the \jTnj) basis, and on the right 
of the equal sign the states are in the \mims) basis. 3? implies "the real part." Given 
'mg = constant' and the symmetry under m; — > —ttli these equations become 



11 
22 



V 



1 1 
22 



1 1 

2 ~ 2 



V 



1 1 
2^2 



i(oi 

3X2 



V 



1 / 1 



oi 

2 



V 



3 \ 2 

3 \ 2 
1\ 1/1 

1 



1-i 

2 

1 

-1 - - 
2 



V 



1 

-1- 
2 



V 



(F.20) 



(F.21) 



Given m,- independence we can take one-half of the sum of these two terms giving 



11 
22 



V 



1/1 
- (0- 
3X2 



1 1 
22 

V 



1 




1 






V 




4)] 


2 


2 





1 

0- 
2 



3 \ 2 



V 



1 / 1 
- - (1 - - 

3 \ 2 

o-i 

2 



V 



1 

1 - - 
2 



2 



1 
3 

-li 

2 



1 

-1 - - 



-1 - - 
2 



[F.22) 



which is equivalent to Eq. ([F.15|) for the specific case of the 2Pi states, as was to be 
shown. 
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